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Abstract 

New lattice data for the Tl u and S~ potentials at short distances are presented. We 
compare perturbation theory to the lower static hybrid potentials and find good agree- 
ment at short distances, once the renormalon ambiguities are accounted for. We use the 
non-perturbatively determined continuum-limit static hybrid and ground state poten- 
tials at short distances to determine the gluelump energies. The result is consistent with 
an estimate obtained from the gluelump data at finite lattice spacings. For the lightest 
gluelump, we obtain Ag s (^/ = 2.5r ( T 1 ) = [2.25±0.10(latt.) ±0.21(th.) ±0.08(A M g)] r 1 
in the quenched approximation with r^ 1 ~ 400 MeV. We show that, to quote sensible 
numbers for the absolute values of the gluelump energies, it is necessary to handle the 
singularities of the singlet and octet potentials in the Borel plane. We propose to sub- 
tract the renormalons of the short-distance matching coefficients, the potentials in this 
case. For the singlet potential the leading renormalon is already known and related 
to that of the pole mass, for the octet potential a new renormalon appears, which we 
approximately evaluate. We also apply our methods to heavy light mesons in the static 
limit and from the lattice simulations available in the literature we obtain the quenched 
result A R£ >/ = 2.5 r^ 1 ) = [1.17 ± 0.08(latt.) ± 0.13(th.) ± 0.09(A M g)] Tq 1 . We calcu- 
late ^bUs( m bMS^ an d a PPly our methods to gluinonia whose dynamics are governed 
by the singlet potential between adjoint sources. We can exclude non-standard linear 
short-distance contributions to the static potentials, with good accuracy. 
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1 Introduction 



In recent years, we have witnessed growing interest in the physics of gluelumps and static 
hybrid potentials. In many cases this has been driven by increasingly reliable lattice sim- 
ulations of their properties [TJ El HI El EH Hj • These results expose models of low energy 
QCD to stringent tests and therefore enhance our understanding of the underlying dynam- 
ics. The short distance physics of the static hybrid potentials is of particular importance. In 
this region, hybrids and gluelumps are intimately related and well suited to investigate the 
interplay between perturbative and non-perturbative physics. At short distances, r, one is 
faced with widely separated scales: l/r> Aqcd- In such situations, effective field theories 
(EFTs) are particularly useful since they enable the physics associated with the different 
scales to be factorized in a very efficient and model independent way. One EFT designed 
to deal with the kinematical case of interest to us corresponds to potential non-relativistic 
QCD (pNRQCD) [S], in the static limit 0. 

In Ref. [HI the gluelumps and the short distance regime of the static hybrids were studied 
within this EFT framework and general features identified. Some results known from the 
past jTUl El HI1 U21 UHj were recovered within a unified framework and in some cases extended. 

One can go beyond this analysis and use lattice data plus the knowledge of the (per- 
turbative) octet potential to obtain numerical values for gluelump masses in a particular 
scheme. However, analogously to the situation with the static singlet potential, the conver- 
gence of the perturbative series of the octet potential does not appear very promising. This 
is a general problem when different scales are factorized, and in particular perturbative from 
non-perturbative ones. The bad convergence is also related to the problem of factorizing 
non-perturbative quantities, without defining their perturbative counterparts [H], and is 
usually believed to be due to the existence of singularities in the Borel transform of the per- 
turbative quantity. These singularities appear to be due to scales of order e~ n x-(the typical 
scale of the perturbative quantity) in an n-loop calculation. In Ref. [T3] one of the present 
authors proposed that, since these singularities are related to energy scales much lower than 
the ones that are supposedly included in the perturbative object, they should be subtracted 
from it and introduced in the matrix elements of the effective theory. This programme has 
been worked out for the pole mass and the static singlet potential [T31 [TE]. Here we apply 
the same approach to the static octet potential. This will allow us to determine absolute 
values for the gluelump masses from the spectrum of the static hybrids, as well as to study 
up to which scale one can use perturbation theory to describe hybrid potentials. 

This paper is organized as follows. In Sec. El we will work out the role of gluelumps in 
pNRQCD, and how gluelumps and hybrid potentials are interrelated. In Sec. El we will then 
sketch how our lattice data have been obtained, before discussing and classifying renormalons 
and power corrections in the continuum MS scheme as well as in a lattice scheme in Sec.|U In 
the same section we will also generalise the renormalon subtracted (RS) scheme of Ref. ^H] 
to the case of the octet potential and discuss the scale dependence. In Sec. E]we will obtain 
the gluelump masses in the RS scheme and relate these results to the lattice scheme. We 
will compare to previous literature and predict the gluelump spectrum. In Sec. |H1 we will 
determine the binding energy of static-light mesons as well as the bottom mass, before we 
discuss generalisations to and relations with adjoint potentials, gluinonium and other objects 
with relevance to short-distance QCD in Sec. 
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2 Hybrid potentials and gluelumps 



We discuss the relationship between hybrid potentials and gluelumps at short distances. 
First we consider the EFT picture, before we discuss the symmetries that are relevant in the 
non-perturbative case. Finally we compare these expectations to lattice data. 

2.1 pNRQCD and gluelumps 

The pNRQCD Lagrangian at leading order in 1/m and in the multipole expansion reads [SI E|, 



^pNrqcd = / d\d 3 RTr js* (ido ~V S )S + Ot (iD - Vo) j - J d 3 R^F* u F^ a + 0(r). (1 



All the gauge fields in Eq. JTJ) are evaluated in R and t, in particular F^ ua = F^ a (R, t) 
and iD O = id O — g[A (R,t),O]. The singlet and octet potentials Vi, i = s,o are to be 
regarded as matching coefficients, which depend on the scale v us separating soft gluons from 
ultrasoft ones. In the static limit "soft" energies are of 0(1/ r) and "ultrasoft" energies are 
of 0(a s /r). Notice that the hard scale, m, plays no role in this limit. The only assumption 
made so far concerns the size of r, i.e. 1/r ^> Aqcd, such that the potentials can be computed 
in perturbation theory. Also note that throughout this paper we will adopt a Minkowski 
space-time notation. 

The spectrum of the singlet state reads, 

E s (r) = 2m os + V s (r) + 0(r 2 ), (2) 

where mos denotes an on-shell (OS) mass. One would normally apply pNRQCD to quarkonia 
and in this case mos represents the heavy quark pole mass. For the static hybrids, the 
spectrum reads 

E H (r) = 2m s + V (r)+A° s + 0(r 2 ), (3) 

where 

A° s = lim i A ln<tf»(T/2)#* (T/2, -T/2)H b (-T/2)) , (4) 



0(T/2,-T/2) = 0(T/2,R,-T/2,R) = Pexp|-^ J^dtA (R,t)\ , (5) 



and H represents some gluonic field, for examples see Table |U in Sec. 15.31 

Eq. (JHJ) allows us to relate the energies of the static hybrids Eh to the energies of the 
gluelumps, 

A H s = {E H (r)-E s (r)}-{V (r)-V s (r)]+0(r 2 ). (6) 

This equation encapsulates one of the central ideas of this paper. The combination Eh ~ E s 
is renormalon-free in perturbation theory [up to possible 0(r 2 ) effects], and can be calculated 
unambiguously non-perturbatively: the ultraviolet (UV) renormalons related to the infrared 
(IR) renormalons of twice the pole mass cancel each other. However, Ah contains an UV 
renormalon that corresponds to the leading IR renormalon of V Q . 

The shapes (of some) of the Eh(t) have been computed on the lattice for instance in 
Refs. [U El El IH El E] • On the other hand the values of (some) Ah have also been computed 
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within a variety of models as well as in lattice simulations [T21EI- Consistency would require 
that the values of A# obtained from Eh — E s and the values of Ah directly obtained from 
gluelump computations should agree. This will be checked in Sec. 15.21 below. 

Gluelump states are created by a static source in the octet (adjoint) representation at- 
tached to some gluonic content (if) such that the state becomes a singlet under gauge 
transformations. This is what would happen to heavy gluinos in the static approximation. 
Hence sometimes gluelumps are also referred to as gluinoballs or glueballinos in the litera- 
ture. Without further information, their energy is only fixed up to a global constant. Only 
the energy splittings between different states have well defined continuum limits in lattice 
simulations. In lattice regularization at a lattice spacing a the normalization ambiguity is 
reflected in a linear divergence oc a -1 while in dimensional regularization one encounters an 
UV renormalon. In the HQET picture of a heavy-light meson one faces a similar problem. 
In this situation one also has a static source (in the fundamental representation in this case), 
which has to be attached to some light-quark (and gluonic) content to become a singlet 
under gauge transformations. The binding energy A is again only defined up to a global 
constant [TH] and only its sum with the pole mass is unambiguous: 



We will investigate this situation in Sec. El below. 

2.2 Symmetries of hybrid potentials and gluelumps 

The spectrum of open QCD string states can be completely classified by the quantum num- 
bers associated with the underlying symmetry group, up to radial excitations. In this case, 
these are the distance between the endpoints, the gauge group representation under which 
these endpoints transform (in what follows we consider the fundamental representation), and 
the symmetry group of cylindrical rotations with reflections D^. The irreducible represen- 
tations of the latter group are conventionally labelled by the spin along the axis, A, where 
E, II, A refer to A = 0, 1, 2, respectively, with a subscript r\ = g for gerade (even) PC = + or 
7] = u for ungerade (odd) PC = — transformation properties. All A > 1 representations are 
two-dimensional. The one-dimensional E representations have, in addition to the rj quantum 
number, a a v parity with respect to reflections on a plane that includes the two endpoints. 
This is reflected in an additional ± superscript. The state associated with the static singlet 
potential transforms according to the representation E+ while the two lowest lying hybrid 
potentials are within the Ii u and E~ representations, respectively. 

In contrast, point-like QCD states are characterised by the J PC of the usual 0(3) <E> C 
rotation group as well as by the gauge group representation of the source. In the pure gauge 
sector, gauge invariance requires this representation to have vanishing triality, such that the 
source can be screened to a singlet by the glue. States created by operators in the singlet 
representation are known as glueballs, octet states as gluelumps. In contrast to gluelump 
states, where the octet source propagates through the gluonic background, the normalization 
of glueball states with respect to the vacuum energy is unambiguous. 

Since C 0(3) ®C, in the limit r — > certain hybrid levels must become degenerate. 
For instance, in this limit, the E~ state corresponds to a J PC = l + ~ state with J z = 




-t-os 




(7) 
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while the Tl u doublet corresponds to its J z = ±1 partners. The gauge transformation 
property of the hybrid potential creation operator will also change in this limit, 3 (g> 3* = 

1 © 8, such that hybrids will either approach gluelumps [cf. Eq. (j2J)] or glueballs, in an 
appropriate normalization. In the case of glueballs the correct normalization can be obtained 
by considering the difference EH{r) — E s (r) from which the pole mass cancels. We will discuss 
the situation with respect to gluelumps in detail in Sec. |U below. 

In perturbation theory, the ground state potential corresponds to the singlet potential 
while hybrid potentials will have the perturbative expansion of the octet potential. 

Recently, Philipsen [TH] suggested to non-perturbatively generalise the octet potential, 
employing a definition that resembles the perturbative one, after gauge fixing to Laplacian 
Coulomb gauge. He proved that this construction is equivalent to a gauge invariant cor- 
relation function whose eigenvalues will resemble masses of physical states. In the limit 
r — > the suggested operator will be an adjoint temporal Schwinger line, dressed with a 
non-local but symmetric gluon cloud, with the J PC quantum numbers of the vacuum. A sim- 
ilar construction is mentioned in paragraph 2 of Sec. El below, as a possible non-perturbative 
normalization point for gluelump energies. The static "octet" potential suggested in Ref. [TJ] 
will have the E+ symmetry and, up to a different non-perturbative off-set, the same pertur- 
bative expansion and power term/renormalon structure as the hybrid potentials discussed 
below. Due to the nature of its creation operator which is non-local, even in the r = limit, 
at present it is not obvious to us how this non-perturbative state can be interpreted in terms 
of the local states we are considering in this paper, certainly an open question that should 
be addressed in the future. 

2.3 Hybrid and gluelump mass splittings 

We would like to establish if lattice data on hybrid potentials reproduces the degeneracies 
expected from the above discussion in the short distance region. In the limit r — > 0, any 
given A > 1 hybrid potential can be subduced from any J state with J > A and PC = + 
for i] = g or PC = — for rj = u representations. For instance the U u is embedded in 
2 + ~, 2 h , • • •. The situation is somewhat different for A = states, which have the 
additional a v parity: the E+ representation can be obtained from ++ , 1 , 2 ++ , • ■ •, E^ from 
, 1 ++ , • • •, E+ from +_ , • • • and E~ from 0~ + , • • •. We list all combinations of 
interest to us in Tabled The ordering of low lying gluelumps has been established in Ref. [T2~] 
and reads with increasing mass: 1 , 2 , 2 H , 3 + ~, ++ , 4 , with a 3 state in 
the region of the 4 and The 2 H and 3 +_ as well as the 4 and 1~ + states are 

degenerate within present statistical uncertainties 1 . The continuum limit gluelump masses 
are displayed as circles at the left of Fig. ^ where we have added the (arbitrary) overall 
constant 2.26/r to the gluelump splittings to match the hybrid potentials. The similarity 
of this value to our estimate of the gluelump energy in Sec. 15. II below is purely accidental. 

Juge, Kuti and Morningstar 1J have, for the first time, comprehensively determined the 
spectrum of hybrid potentials. We convert their data, computed at their smallest lattice 
spacing a a pa 0.2 fm, into units of tq pa 0.5 fm [20]. Since the results have been obtained 

1 The splittings between all states with respect to the l" 1 ground state have been extrapolated to the 
continuum limit in Ref. |12j and we add our own extrapolations for the 4 and 1 ^ states to these, based 
on the tables of this reference. 
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point particle J FC 


open string A° v 


1+- 




1~ 




2— 




2+- 




3+- 




0++ 




4~ 




1-+ 


y+r ttw 



Table 1: Expected degeneracies of hybrid potentials at short distance, based on the level 
ordering of the gluelump spectrum. Note that if the 3 + ~ gluelump turned out to be lighter 
than the 2 + ~ then the E~ , IT4, A u , $ u potentials would approach the 3 +_ state while the 
£+,H"", A^ potentials would approach the 2 + ~ instead. 

with an improved action and on anisotropic lattices with a T w a CT /4, one might expect lattice 
artifacts to be small 2 , at least for the lower lying potentials. Hence we compare these data, 
normalized to i? s +(ro), with the continuum expectations of the gluelumps The full 

lines are cubic splines to guide the eye while the dashed lines indicate the gluelumps towards 
which we would expect the respective potentials to converge. 

The first 7 hybrid potentials are compatible with the degeneracies suggested by Table [T] 
The next state is trickier since it is not clear whether 2 + ~ or 3 + ~ is lighter. In the figure we 
depict the case for a light 2 + ~. This would mean that (E+, IT4, A u ) approach the 2 +_ while 
(E~', IT", A^, $ u ) approach the 3 + ~. Note that of the latter four potentials only data for JJ" 
and Q u are available. Also note that the continuum states 174, and $ u are all obtained 
from the same E u lattice representation. For the purpose of the figure we make an arbitrary 
choice to distribute the former three states among the E' u , E'^ and lattice potentials. To 
firmly establish their ordering one would have to investigate radial excitations in additional 
lattice hybrid channels and/or clarify the gluelump spectrum in more detail. Should the 
2 + ~ and 3 H hybrid levels be inverted then (£„', IT4, A u , $ u ) will converge to the 3 + ~ while 
(£+, IT", A' u ) will approach the 2 + ~. We note that the ordering of the hybrid potentials, with 
a low £+, makes the first interpretation more suggestive. 

Finally the £+" potential seems to head towards the ++ gluelump but suddenly turns 
downward, approaching the (lighter) sum of ground state potential and scalar glueball |2T| l2^] 
instead. The latter type of decay will eventually happen for all lattice potentials but only 

2 On the lattice the relevant symmetry group is D^h rather than Doo/i (see e.g. Ref. [23j)- In the continuum 
limit the potentials will correspond to £+, the A2 V potentials to S~ and the E v potentials to 11^, where 
rj = u, g. The radial excitations could in principal correspond to higher spin potentials and in fact one of the 
three observed excitations of E u will correspond to the <&„ ground state. In all other cases, associating the 
lowest possible continuum spin to a given lattice potential seems to agree with the ordering suggested by the 
gluelump spectrum (as well as in the large distance string limit pQ). Bi^ and B277 both correspond to A v . In 
either case (as well as for A^), at the short distances displayed in the figure, the two lattice representations 
agree with each other, supporting the view that violations of rotational symmetry are small. In this case we 
only display the lattice representations with better statistical accuracy, i.e. the b['^s. 
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0.5 1 1.5 2 2.5 

r/r 

Figure 1: Different hybrid potentials J^j at a lattice spacing a a ~ 0.2 fm ~ 0.4 ro, where ro ~ 
0.5 fm, in comparison with the gluelump spectrum, extrapolated to the continuum limit \1 6 2\] 
(circles, left-most data points). The gluelump spectrum has been shifted by an arbitrary 
constant to adjust the l + ~ state with the IT U and S~ potentials at short distance. In addition, 
we include the sum of the ground state ) potential and the scalar glueball mass m ++ JUTl 
22]. The lines are drawn to guide the eye. 

at extremely short distances. We also remark that all potentials will diverge as r — ► 0. This 
does not affect our comparison with the gluelump results, since we have normalized them to 
the n u /E~ potentials at the shortest distance available. (The gluelump values are plotted 
at r = to simplify the figure.) 

On a qualitative level the short- distance data are very consistent with the expected 
degeneracies. From the figure we see that at r w 2 r ~ 1 fm the spectrum of hybrid 
potentials displays the equi-distant band structure one would qualitatively expect from a 
string picture. Clearly this region, as well as the cross-over region to the short-distance 
behaviour r < r < 2r , cannot be expected to be within the perturbative domain: at best 
one can possibly imagine perturbation theory to be valid for the left-most 2 data points. With 
the exception of the U u , U' u and $ u potentials there are also no clear signs for the onset of 
the short distance 1/r behaviour with a positive coefficient as expected from perturbation 
theory. Furthermore, most of the gaps within multiplets of hybrid potentials, that are to 
leading order indicative of the size of the non-perturbative r 2 term, are still quite significant, 
even at r = 0.4 ro; for instance the difference between the S~ and IT U potentials at this 
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smallest distance is about 0.28 r^ 1 « 110 MeV. 



2.4 The difference between the II u and E w hybrids 

From the above considerations it is clear that for a more quantitative study we need lattice 
data at shorter distances. In this paper we have obtained these for the lowest two gluonic 
excitations, U u and (see Sec. OJ). We display their differences in the continuum limit 
in Fig. 121 We see how these approach zero at small r, as expected from the short distance 
expansion. pNRQCD predicts that the next effects should be of 0(r 2 ) (and renormalon-free) . 
In fact, we can fit the lattice data rather well with a AE U v + = A u v -r 2 ansatz for short 
distances, with slope (see Fig. |2J), 

A nu _ 2 -=0.92l°ir 3 , (8) 

where the error is purely statistical (lattice). This fit has been done using points r < 0.5 r$. 
By increasing the fit range to r < 0.8 r the following result is obtained, 

K u -v- = (0-83 ± 0.29) r" 3 , (9) 

indicating stability of the result Eq. (jHJ) above. 

In order to estimate systematic errors one can add a quartic term: br A (only even pow- 
ers of r appear in the multipole expansion of this quantity). If the result is stable, our 
determination of ^4n„-s~ should not change much. Actually this is what happens. If we 
fit up to r < 0.5 r , we obtain the central value y4n„-s u ^o = 0-93 with a very small quartic 
coefficient, br$ = —0.05. If we increase the range to r < 0.8 r , we obtain the same central 
value, ^4n u -£- r o = 0.93, but with a slightly bigger quartic term, br$ = —0.18. Introducing 
the quartic term enhances the stability of A u _ s - under variations of the fit range. From 
this discussion we conclude that the systematic error is negligible, in comparison to the error 
displayed in our result Eq. (jSJ). 

We remark that within the framework of static pNRQCD and to second order in the 
multipole expansion, one can relate the slope ^4 nu _ s - to gluonic correlators of QCD. 



3 Lattice determination of hybrid potentials 

We extract the hybrid potentials in two sets of simulations, using the Wilson gauge ac- 
tion on an isotropic lattice with volume 24 3 x 48 at (3 = 6.2 (a ~ 0.14 r ) as well as 
on three anisotropic lattices with spatial lattice spacings a a ~ 0.33, 0.23, 0.16 r , respec- 
tively, with anisotropy a a « 4a T . The former result has been obtained in the context 
of the study of Ref. |I] (and has been published in Ref. [3 ) while the simulation pa- 
rameters, statistics and smearing of the latter runs are identical to those of Ref. j24j : 
(/^>£o) = (5.8,3.1), (6.0,3.2), (6.2,3.25). The isotropic data are used as a consistency check 
and in Sec. 14.41 below, while we extrapolate the data obtained on the anisotropic lattices to 
the continuum limit. 

Some time was spent on improving the shape of the hybrid creation operators to optimize 
the overlap with the ground state [3] . The Il n potential has been determined on-axis as well 
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0.2 0.4 0.6 0.8 1 1.2 1.4 

r/r 

Figure 2: Splitting between the £~ and the U u potentials, extrapolated to the continuum 
limit, and the comparison with a quadratic fit to the r < 0.5 r data points (tq ~ 0.4 GeV). 
The big circles correspond to the data of Juge et al. JJ], obtained at finite lattice spacing 
d(j 0.39 r . The errors in this case are smaller than the symbols. 

as along a plane-diagonal, r/a CT oc (1, 1,0), while the potential has only been obtained 
on-axis. Typically we achieved ground state overlaps of around 65 % for both potentials 
at (3 = 5.8 and between 85 % and 90 % at the larger two (3 values. Typical fit ranges for 
one-exponential fits to correlation functions for the U U (Y>~) potential were 8 < t/a T < 18, 
(9 < t/a T < 14) at (3 = 5.8, 9 < t/a T < 24 (11 < t/a T < 21) at /3 = 6.0 and 13 < t/a T < 30 
(15 <t/a T < 25) at (3 — 6.2. For all further details of the analysis we refer to Ref. [21] where 
potentials between sources in non-fundamental representations of SU(3) were extracted using 
exactly the same methods. 

Subsequently, the potentials as well as differences between potentials have been extrapo- 
lated to the continuum limit. As one such example we display the difference between the Il n 
and the singlet potential in Fig. El for distances r < tq. In this extrapolation we somewhat 
deviate from Ref. [21]: we follow Ref. [22] in removing the lattice artifacts to leading order 
in a s , by plotting the data as a function of the inverse lattice Coulomb propagator, 

r i 



(10) 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

f/r 

Figure 3: Continuum limit extrapolation of the difference between the U u and the poten- 
tials vs. r = r[l + 0(a 2 /r 2 )] as described in the text [Eqs. Hu}) and < f77)) /. The Juge et al. 
data are from Ref. 



rather than of r. The lattice Coulomb propagator for the Wilson gauge action is given by. 

3 Q cos(QR) 



■ 1 - 






= 4tt f 


R_ 


L J- 



-(2tt) 3 4 Ei sin (f) 



(11) 



and agrees with the continuum l/i?-function up to 0(a 2 /r 2 ) lattice artifacts. R = r/a 
denotes an integer valued three-vector and the Qi = qia a are dimensionless. For the H u 
potential this procedure removes violations of rotational symmetry within the statistical 
errors and brings the plane- diagonal points in-line with the on-axis data. Unfortunately, we 
cannot perform a similar internal test for the S~ potential which we only determined for 
on-axis separations. 

The next step involved fitting differences between hybrid potentials and AEh = 
Eh — E^+ , for r > 2a to the phenomenological interpolation, 



AfaH^.+ ^^M +o.r', (12) 

r 

with parameters q. We then extrapolated these interpolating curves to the continuum 
limit, assuming the leading order a 2 dependence. This was done separately for different 
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pairs of two lattice spacings. The central value of the extrapolation is given by the result 
obtained from the a a ~ 0.33 r and a a ~ 0.16 r data sets. The error is estimated by the 
squared sum of the statistical error of the fine lattice data set and the difference between the 
above extrapolation and an extrapolation obtained from the a a ~ 0.23 ro and a a « 0.16ro 
data sets. With decreasing r the interpolating fits become less well constrained and hence 
the latter systematic uncertainty increases. The resulting error band is depicted in Fig. El 
Reassuringly, the a a « 0.16 r data are already in agreement with the continuum limit 
and the a a « 0.23 r data agree within errors: the fine lattice data set effectively already 
corresponds to the continuum limit. The more precise isotropic reference data (a ~ 0.14 ro) 
are also close to the continuum limit. We also notice that the first three data points of the 
coarse lattice data by Juge et al. lj (<V ~ 0.39 r ) are compatible with our extrapolation. 
The same observations hold true for the E~ potentials. 

Rather than representing the continuum limit extrapolated potentials by error bands, in 
the remaining parts of this paper we add the difference between (finite a) interpolation and 
(continuum limit) extrapolation to the fine lattice data points and increase their errors by 
the systematic uncertainty involved in the extrapolation. 

4 Static octet potential 

We will discuss the octet potential in the OS (= "pole mass") scheme, compute the normal- 
ization constant of the renormalon and generalise the RS renormalon subtracted scheme 
to this case. We will also discuss the structure of power divergences on the lattice and the 
analogous lattice scheme. Finally we discuss the running of the gluelump mass from one 
scale to another. 

4.1 OS scheme for the octet potential 

The octet potential in the case l/r> A QCD can be computed order by order in perturbation 
theory. Nevertheless, it is not an IR safe object j2H]- Its perturbative expansion reads, 



where we have made explicit its dependence on the IR cutoff v us and a s = a s (^), where we 
define 



In what follows we will always identify a s with oc^s. The first two coefficients V 0i q, V _\ are 
known, as well as the leading-log terms of V $ (for the renormalizat ion-group improved 
expression see Ref. |2~7j). Note, however, that these leading logs are not associated to the 
first IR renormalon. For V 0j 2 there exists a preliminary computation j2H] , 



oo 



V {r; v us ) ~ V o,nOt n s 



(13) 



n=0 




V , 2 



1 



K,2 + 5V ,2 , 



(14) 



i\T c 2-l 
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5V , 



1 :21Cl- 



(15) 



2N C {Air) 2 ' 

which we will use in what follows. V s ^ has been computed in Ref. For V Q $, we will use 
the renormalon-based estimate that we obtain in Sec. 14.21 below (Table EJ). 
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Figure 4: roV (r) (the octet potential in the OS scheme) at tree level (dashed lines), one-loop 
(dashed- dotted lines), two loops (dotted lines) and three loops (estimate) plus the leading 
single ultrasoft log (solid lines). Fig. a) corresponds to the scale v = z/j [cf. Eq. Mb}) ] and 
Fig. b) to v = 1/r. In both cases, v us = 2.5r (j 1 . Only the solid curves depend on this choice. 

Studying the convergence of perturbation theory of the octet potential in the OS scheme, 
conclusions similar to those in Ref. ^H] are obtained. The poor convergence is demonstrated 
in Fig. HI where we try two choices of the scale v. In part (a) we use v — Ui where, 



j-i 



6.604 r^ 1 



2.6 GeV. 



(16) 



corresponds to the shortest distance r' for which the continuum limit extrapolated lattice 
potentials are available. In part (b) we vary v = 1/r. Obviously the curves depicted in the 
two parts of the figure agree with each other at r = r' ~ 0.15 r$. Note the difference in the 
vertical scale. 

4.2 Static octet potential normalization constant 

We define the Borel transform of the octet potential as follows, 



V = I dte-^BlV^t) 



B[V }(t) = J2V ,n-. 



n=0 



(17) 



The behaviour of the perturbative expansion Eq. (|T3|) at large orders is dictated by the closest 
singularity to the origin of its Borel transform, which happens to be located at t = 2tt/13q. 
This singularity has two sources: one is an UV renormalon which cancels with the renormalon 
of twice the pole mass, the other is an IR renormalon that cancels with the UV renormalon 
of the gluelump energy. This result follows from the structure of the effective theory and the 
consequent factorization of the different scales in Eq. (jH)). Being more precise, the behaviour 
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of the Borel transform of the static octet potential near the closest singularity to the origin 
[u = 1/2 where we define u = j3 t/(4:Tr)] reads, 

B[V ](t(u)) = N Vo" ^_ 2u y +b (l + ci(l - 2u) + c 2 (l - 2uf + ■■■) + (analytic term) , 

(18) 

where by analytic term, we mean a function that is analytic up to the next IR renormalon 
at u — 3/2. This dictates the behaviour of the perturbative expansion at large orders to be, 

v W rtn + 1 + 6) / 6 6(6 - 1) \ 

K,„ = Wv.^ s J r(1 + >) ( 1 +(^j*+ ( ,, + t)( ,, + t _ 1 ) '> + -J- d») 

The structure of the renormalon is equal to the singlet one. This is due to the fact that the 
number of octet fields is conserved at leading order in the multipole expansion and that the 
mass (potential) does not renormalize at this order. Therefore the values of the coefficients 
6, ci, C2, • • • above are the same as for the case of the static potential and the pole mass and 
can be found in Refs. (SHI EH • We display them here for ease of reference: 



2 



(20) 



46/3 3 U 



Ci = TT5al^-A) » (21) 



and 



_ 1 (3j + 4/3 3 A& - 2(3 (3l(3 2 + (3l{-2(3\ + /g) - 2/3 4 /3 3 
c 2 - 171 TT ^TTTis ■ 



6(6 - 1) 32/3, 



The only difference with respect to the static singlet potential is the value of Ny . The 
cancellation of the renormalon in Eq. (J3J) requires, 

2N m + N Vo + N A = 0, (23) 

where TVa is the normalization constant of the renormalon of the gluelump mass (B[A] 
reads the same as Eq. (|18j) . with the replacement Ny o i— > N\). Therefore, unlike in the 
static singlet potential case, we cannot fix Ny o from the knowledge of N m alone. Yet we will 
(approximately) determine Ny o from low orders in perturbation theory of the octet potential. 
Note also that N\ is independent of H, the specific gluonic content of the gluelump, since it 
only depends on the high energy behaviour, which is universal. To leading non-trivial order 
one obtains, N Vo = C A /2 - C f , N A = C A /2. 

In analogy to Refs. |321 HH1 UHj we define the new function, 

oo 

DvM = ^4:V = (1-2 M ) W B[K (0) ](^)) (24) 

n=0 

= N Vo v (l + cj(l - 2u) + c 2 (l - 2m) 2 + ••■) + (1 - 2n) 1+b (analytic term) , 

and try to approximately determine Ny a by using the first three coefficients of this series. In 
analogy to Refs. [EHE)], we fix v = 1/r and obtain [up to 0(u 3 )\ u= i/2], 

N Vo = 0.166667 - 0.0624292 + 0.00976333 = 0.114001 . (25) 
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The convergence is rather good and, moreover, we have a sign alternating series. In fact, the 
scale dependence is becoming milder when we go to higher orders (see Fig. |5J). Note that 
if the two loop coefficient V^ 2 had been equal to that of the singlet case [2H| (with colour 
factor Cf i — > Ca/2 — C/), we would have obtained Ny = 0.146542. 
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Figure 5: x = vr dependence of Ny a for rif = at LO (dashed- dotted line), NLO (dotted 
line) and NNLO (dashed line). 

We can now compute estimates for V 0iTl by using Eq. (fTTJj) . These, as well as estimates for 
V s>n , are displayed in Table El for n/ = 0. We can see that the exact results are reasonably 
well reproduced. Hence we feel confident that we are near the asymptotic regime dominated 
by the first IR renormalon and that for higher n our predictions will accurately approximate 
the exact results. 



V = rV 

' o,n 1 v o,n 


V ,o 


Vo,l 


V 0t2 


V 0>3 


Vol 


exact (rif = 0) 
Eq. (HHD {n f = 0) 


0.166667 
0.110552 


0.305472 
0.244266 


1.27419 
1.14193 


6.97413 


54.4562 


V = rV 

v s,n 1 ' s,n 


V s>0 


Vsl 


V s , 2 






exact {rif = 0) 
estimate (nf = 0) 


-1.33333 
-1.20643 


-2.44378 
-2.66564 


-11.7893 
-12.4616 


-76.1075 


-594.2718 



Table 2: Values of V a _ n with v = 1/r: exact result (where available) and the estimate using 
Eq. \19\j . We also display the estimates ofV s ^ n with v = 1/r (extracted from Ref. ^B|)■ 

In order to avoid large corrections from terms depending on u us , the predictions should 
be understood with v us = 1/r and later-on one can use the renormalization group equations 
for the static potential |2Zj to keep track of the v us dependence. 
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4.3 RS scheme for the octet potential 

In Sec. 14.11 we have demonstrated the poor convergence of the perturbative expansion of 
the octet potential in the OS scheme. This bad behaviour is usually believed to be due to 
the singularities in the Borel transform of the perturbative expansion. Nevertheless, these 
singularities are fake since they cancel with singularities in the matrix elements. On the 
other hand this lack of convergence of perturbation theory arises because at higher orders in 
perturbation theory smaller and smaller momenta contribute to the short- distance matching 
coefficients of the effective theory. This clashes with the logic of scale separation in the 
EFT formalism. The solution advocated in Ref. ^3] was to subtract this behaviour from 
the matching coefficients. At the practical level this was implemented by subtracting the 
Borel plane singularities of the matching coefficients. In Refs. [To"! ITo] this has been worked 
out for the pole mass and the static singlet potential and we refer to these references for the 
definitions and further details. In particular Eq. (J5J) reads, 

E s (r) = 2m RS (u f ) + V, jBS (r; u f ) + 0(r 2 ) , (26) 

where 

mns(vf) = moa - Smasfa) , (27) 

K,R S (r; vj) = V s (r) + 26m BS {u f ) , (28) 

and (in the above equation we have already used the fact that the renormalon of the singlet 
potential cancels with the one of minus twice the pole mass 



,3 



A ( \ S^at (P°Y n+u T(n+l + b-k) 

8m RS (u f )=^N m v f a s + („,) gc fe r{1 + b _ k) • (29) 

For the static hybrids, the spectrum reads, 

E H (r) = 2m RS (z/ / ) + T/ 0>RS ( r ; u f ) + Af(u f ) + 0(r 2 ) . (30) 

Obviously, we have to define the octet potential and the gluelump mass above. In the RS 
scheme the octet potential reads, 

oo 

VoM"f) = V o- 5V o,KS = E K R n«s n+1 , (31) 



ra=0 



where 



°' RS = S ° f W " s {Uf) £ Ck m+b-k) ■ (32) 

This specifies the gluelump mass which reads, 

A| S (^/) = A ff - 5A RS (u f ), (33) 

3 Actually, throughout this paper we use the RS' scheme as defined in Ref. ^H] instead of the RS scheme, 
since we believe this to have a more physical interpretation. For simplicity of notation we will however refer 
to this modified scheme as "RS scheme" , omitting the "prime" . 
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Figure 6: r V o ^s a t t ree level (dashed lines), one-loop (dashed- dotted lines), two loops (dotted 
lines) and three loops (estimate) plus the leading single ultras oft log (solid lines). For the 
scale of a 8 (u), we set u — (stable behaviour at large distances) or v = l/r (diverging at 
large distances). We kept Uf = 2.5 r^ 1 fixed. 



where 

5A RS (u f )=^N AH , f d?\v t ) Ec fc V(i + & -jfc) • ^ 

Note that the potentials and A^ s depend on Uf which, in the context of pNRQCD, can be 
thought of as a matching scale between ultrasoft and soft physics. In what follows, we will 
set Uf = 2.5 r^ 1 . Results for different values of Uf can be obtained using the running on Uf, 
which is renormalon independent. 

Analogously to the discussion of Ref . jTB] , we can study the convergence of the perturba- 
tive expansion in the RS scheme. In Fig. El we can see that the stability is greatly improved, 
compared to the OS scheme discussed in the previous section. No matter whether we choose 
to work with a s (Vj) or a s (l/r), the expansions converge towards the same curve. In Fig. [7| 
we can also see that they agree with the continuum limit lattice data (we have to subtract 
an unknown constant for this comparison). In this figure the errors of En u (r) — En u (r') for 
r > r' are purely statistical while the (strongly correlated) systematic error of the continuum 
limit extrapolation is only displayed for the first data point [En u {r') — En u (r') = 0], where 
it is largest. 

The price we pay to obtain convergent expansions in a s for the potentials is the introduc- 
tion of power-like terms (proportional to Uf, with logarithmic corrections). This behaviour 
very much resembles that of lattice regularization with a hard cut-off which we discuss below. 
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Figure 7: r [V Ot Yis( r ) ~ K,Rs( r ')] + C at tree level (dashed lines), one-loop (dashed- dotted 
lines), two loops (dotted lines) and three loops (estimate) plus the leading single ultras oft log 
(solid lines) compared with the non-perturbative continuum- limit results for Eu u (r) —En u (r') 
(symbols with error bars). For the scale of a s {v), we set v — 14 = 1/r' (stable behaviour at 
large distances) or v = 1/r (diverging at large distances). A (small) constant C is arbitrarily 
adjusted to show agreement with the lattice data. 



4.4 Lattice scheme for the octet potential 

It is conceptionally illuminating also to consider the situation in lattice regularization. In 
this case, the inverse lattice spacing a -1 results in a hard UV cut-off of the gluon momenta. 
Feynman diagrams are UV finite and EFT matrix elements are manifestly renormalon-free 
as long as they are obtained in non-perturbative numerical simulations. The price paid is the 
existence of power divergences oc a -1 , which cannot be eliminated in the continuum limit. 

The analogy with the previous sections can be made quite evident. In particular, all 
the quantities that we have defined in the OS and RS schemes can also be defined in a 
lattice scheme. There are some differences however. The lattice gluelump A#(a) has a 
power divergence to start with (which can be traded in for a renormalon ambiguity when 
subtracted in perturbation theory). In this sense it is similar to A^f(uf). While formally 
many expressions resemble those of the RS case, a" 1 plays a slightly different role than z/y 
that separates soft from ultrasoft scales since a < r <C i^j 1 . Another difference is that at 
finite lattice spacings the potentials remain finite as r — > 0. In particular, we will see that 
gluelumps are the r — > limits of hybrid potentials (at finite lattice spacing), in perturbation 
theory as well as non-perturbatively. This should not be surprising since the r — ► limit at 
finite lattice spacing corresponds to the situation r « a. This means that the ultraviolet 
cutoff ~ a -1 is much smaller than r _1 and that the dynamical degrees of freedom are only 
the ultrasoft ones. Actually, in this situation, i/f and 1/a play an analogous role. 

Let us illustrate the above by first considering perturbation theory, before discussing the 
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scale separation and how the lattice scheme translates into other schemes, at finite lattice 
spacings as well as in the continuum limit. 

For simplicity we will consider the Wilson discretization of the continuum action. In 
this case the "lattice Coulomb term" [1/R]^ takes the form Eq. (fTTj) . For instance one can 
calculate the finite value, [1/0] £ = 3.17. . .. Using this notation, one finds the lattice results, 



—CfO>La 1 
C A 



K,i(r;a) = 
where the "self energy" is given by, 



1 

R 



[l + 0(a L ))+26rn^ t (a) 



Cf I OtLCL 



L 

-1 



1 

R 



[1 + 0(a L )) + 25ml 



(35) 
(36) 



a <^ s tat(a) 



C f 



-Ol L 



Cf ( 



+ ■ ■ 



+ 



(37) 



Note that unlike in dimensional regularization, by using a hard cut-off, such power diver- 
gencies appear naturally as part of the perturbative expansion. Eqs. (J33j) and (}3l)|) are both 
known to 0(af) and Eq. (as well as the difference V Q — V s ) is also known approximately 
to 0(a^), up to 0(a^a 2 /r 2 ) lattice corrections jHS]- In pure gauge theory with Wilson action, 



(38) 
(39) 
(40) 



the coefficients of the expansion of 5m s L tat read [Ml CHI ESI ES] , 



v 3 



3.1759115..., 
0.21003(5) x 10 3 
20.4(3) x 10 3 . 



= 3/(27r /3) denotes the lattice coupling at a scale a 1 which can be translated into other 

(41) 



schemes such as MS by means of a perturbative computation 



Oil 



a H [a 



bv 



Ali- 



bi - 2bi 



a 2 (a x ) 



(4tt) 



+ 



with [3*6] . 



73.93539066 . . . 
b\ + 1388.1645. 

Let us now consider the singlet case. We have, 



bi 



E% + (t; a) = V s , L {r- a) + A QCD 0(A 2 QCD r 2 ) + 0(A 2 QCB a 2 ) + 0{a 2 /i 



(42) 
(43) 



(44) 



where Aqcd represents a generic non-perturbative scale like r^ 1 . The last two terms account 
for possible non-perturbative lattice artifacts, which vanish as a — > 0. From the quarkonium 
energy E s (r) at r ^> a, we can non-perturbatively obtain the heavy quark mass in a lattice 
scheme, 



m L (a 



1 r 

2 



E s (r)-E^ + ( r] a) +0(a 2 /r 2 ) 



(45) 
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By redefining, 

F s , L (r; a) = V s , L (r; a) - 25m s L tat (a) , (46) 
we can then achieve formal correspondence to Eqs. (|26jl and respectively: 

£ s (r) = 2m L (a) + K, L (r;a) + 0(r 2 ) (47) 
= 2m os +F s , L (r;a)+0(r 2 ), (48) 

where the above two equations are correct up to 0(AQ CD a 2 ) and 0(a 2 /r 2 ) lattice corrections. 
We can relate the heavy quark mass in the lattice scheme to the OS scheme, 

m L (a) = m os - <*mf tat (a) . (49) 

^ m stat( a ) contains the same renormalon as m os , such that Eq. has good convergence 
properties when expanded in terms of a s . m L (a) is proportional to a -1 , with logarithmic, 
as well as 0(a 2 ) lattice corrections. One can convert m^(a) order by order in perturbation 
theory into say Tn^(i/), without renormalon ambiguity. 

In the lattice scheme we also have E^+(0;a) = V s z,(0;a) = 0: the sources are "smeared 
out" on a scale a since the gluon, due to the UV cut-off, cannot resolve structures smaller 
than the lattice spacing. Consequently, the Coulomb term does not diverge as r — > but 
approaches a finite value in units of a. In perturbation theory, in the limit r — > 0, the 
lattice [1/R]^ term exactly cancels with 25m^ tat : the perturbative expansion of Vg^rja), 
Eq. ()35|) above, does not contain the renormalon associated with the pole mass. Non- 
perturbatively, in the limit r — > the Wilson loop becomes a time independent constant, 
such that E^+(0; a) = too. As r > the perturbative V s l acquires a power term. 

Next we consider the hybrid case. We can calculate the gluelump mass in perturbation 
theory 4 , 



a5A L (a) = ~^oi L 



+ ... = c± aL ( Vl + V2 ?*) , . m 



l 2 " V 4tt, 

where v\ and v 2 are the same as for the case of <5m^ at and can be found in Eqs. (|55j) and (J39j) 
above. Note that the 0(a^) term is expected to be different and is not known at present. 
However, A H is related to the difference between V Q and V s , such that any difference with 
respect to the v 3 of Eq. (|3U|) above will be suppressed by a colour factor l/N 2 . 

The tree level expression for V ,l is displayed in Eq. (J3BJ). While the perturbative ex- 
pansion of V Sj l was unaffected by the renormalon of the pole mass, the one of V 0> l contains 
the same renormalon as the expansion of SA^. For r ^> a the renormalon-free combi- 
nation V^x(r;a) — 5A L (a) plays the role of V ^s(r; Vf) in Eq. (|5U|1 . At r = we have, 
V o l(0;o) = 6A L (a) as well as the non-perturbative equality, 



E% u (0;a)=E% (0;a)=A£(a). (51) 



We redefine, 



V 0>L {r; a) = V QtL (r; a) - 8A L {a) - 25m^ at (a) , (52) 



In the context of perturbation theory we do not distinguish between different gluelumps since the mass 
splittings have an entirely non-perturbative origin. 
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Figure 8: Splitting between the lowest two hybrids and the £+ potentials (pentagons and 
squares) as a function ofr/a [see Eq. KTDjl ] at a fixed lattice spacing, a ~ 0.137r 0; in 
comparison to V 0jL (v; a) — V S:L (v; a) = V 0i l(t; a) — V 0i L(r; a) + 5A L (a) at tree level (dashed 
lines, diamonds), one loop (dashed- dotted lines, circles) and two loops (dotted lines, r 3> a 
estimates). The open symbols correspond to the respective gluelumps, non-perturbatively 
(square with pentagon) and in lattice perturbation theory (diamonds and circles). 



to achieve formal correspondence with Eqs. Q and ©: 

E L H (r; a) - E^+(r; a) = A L H {a) + \V , L (r; a) - F s , L (r; a)l + 0(r 2 ) . 



(53) 



Note that Ejj(r; a) = Eh(t) — 2m^(a) + 0(a 2 /r 2 ), in analogy to Eq. (|45p. The combination, 



V 0)L (r; a) - V ' BjL (r; a) = ^-a L a~ 



'( 


■ 1 ■ 




V 


J 




R_ 


L 


.0. 





+ 0(a 2 



(54) 



vanishes for r = and is renormalon-free. The same holds true for Ejj(0; a) — E^ + (0;a) 
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Aff(a) = 0: Eq. is not only valid for r > a but also for 5 r = 0. We have, 

A L H (a) = A° s + 5A L (a) . 



(55) 



Note that the above equation is only correct up to non-perturbative 0(Aq CD ci 2 ) contributions 
to A^r . Again A^ is renormalon-free but has a power divergence. By subtracting SAi(a) 
order by order in perturbation theory one can obtain an on shell A^ s , but at the price of a 
renormalon ambiguity Note the similarity between the above equation and Eq. (pI3j) . 
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Figure 9: Splitting between the lowest two hybrids and the potentials (pentagons and 
squares) as a function of r/a [see Eq. KWi) } at a fixed lattice spacing, a ~ 0.137r , in 
comparison with K oi (r;a) — V 0i L(r;a) + C at tree level (dashed lines, diamonds), one loop 
(dashed- dotted lines, circles) and two loops (dotted lines, r 3> a estimates). The vertical 
normalization C has been adjusted to produce agreement at r/a = a/2. The open symbols 
correspond to the respective gluelumps, non-perturbatively (square with pentagon) and in 
lattice perturbation theory (diamond and circles). 

5 Based on the results of Sec. 15.21 below as well as of Ref. we know that the l" 1 glueball will become 
lighter than the gluelump Aj^(a) around a < r c w tq/7, when using the Wilson action. In fact we discussed 
a similar situation in Sec. 12.31 above, for the £+' potential. This limit is not yet relevant for the U u and 
E~ potentials at the lattice spacings covered in this paper. In the case a < r c , Eq. (|53|) will still apply for 
t~ 3> 3> Aqcd, however, Eq. Q51J1 will become modified; it would apply to the first radial excitations 
in the hybrid channels rather than to the ground states, until finally around a ~ f*o/12 a continuum of 
two-glueball states is encountered. 
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In Fig. |H1 we compare non-perturbative data on the splitting between hybrid potentials 
with respect to the ground state potential with the perturbative expectation. The data have 
been obtained by us on an isotropic lattice at (3 — 6.2 with lattice spacing a ~ 0.137 r . 
Both gaps, Eh — (squares) and Ez,- — E^+ (pentagons) are plotted as a function of 

r/a [see Eq. (|10j) ]. The differences are indicative of the size of the expected non-perturbative 
0(r 2 ) contributions. We compare the non-perturbative data to the perturbative expecta- 
tion for V 0) l{y'i a) — K,l(f;gO- The latter perturbation theory will suffer from the same 
renormalon ambiguity as 0Al(ci) and the difference between perturbation theory and non- 
perturbative data corresponds to A^ s . The left-most points (open symbols) correspond to 
the A B gluelump, plotted at r/a = [1/0]l ~ 0.315. 

The evaluation was done both in terms of a s (a~ 1 ) and a s (uf) where Vf = 2.5 r^ 1 « 
0.34 a -1 . To simplify the figure we disregard the uncertainty in the determination of A^g = 
0.602(48) r^ 1 [37]. At LO and NLO lattice perturbation theory results are available 13*31 
(diamonds and squares). Since everything is plotted as a function of r/a = [1/R]^ 1 all 
diamonds lie exactly on top of the dashed continuous r>a curves while at small distances 
there are differences between the dashed-dotted NLO curves and the exact NLO results 
(circles). In addition we plot the r 3> a limits to NNLO (dotted curves). The shapes of the 
perturbative curves remain qualitatively stable while the normalization is not convergent as 
the order of the expansion is increased and is also strongly affected by the scale of a s {v). 
This behaviour reflects the presence of the renormalon of A§ s , quite similar to what we can 
see in Fig. @Jl 

By comparing with the renormalon- free right hand side (rhs) of Eq. ()53j) a better conver- 
gence can be achieved. However, such a comparison is only possible up to 0(a 2 ) as we do 
not exactly know the 0(al) contribution to the counterterm 0Al(ci) in the lattice scheme. 
Instead we choose to demonstrate the quality of the perturbative expansion in Fig. |U] by 
adding global normalization constants to all curves in such a way that agreement is pro- 
duced at r/a = \/2. (We shall return to the question of renormalon cancellation in A^-(a) 
in Sec. 15.21 below.) Indeed the differences between NNLO and NLO are smaller than those 
between NLO and LO. Moreover, at higher orders the scale dependence is reduced. The 
v = a^ 1 curves seem to describe the data better at small r while the v — Vf curves seem 
to work better at larger r. Up to distances as big as f = r ~ 7.3 a the perturbative curves 
seem to have an accuracy better than the non-perturbative uncertainties, estimated by the 
difference _E s -(r) — En u {r). 

We mentioned above that while formally the lattice spacing a~ l appears in the same 
places in the lattice scheme as the scale Uf did in the RS scheme of Sec. 14.31 these two scales 
should not be confused with each other as a -1 > r _1 > i/f > Aqcd- Conceptionally we have 
been discussing the situation in which the potentials are evaluated in perturbation theory 
at scales v > i/f while A^- s is an ultrasoft matrix elements, associated to physics at scales 
smaller than i/f. The lattice encapsulates the same physical picture. For instance, to each 
finite order in perturbation theory 6 , V s / 0> l{y; a) r ^5 > 25m^ at (a) and V^l(0, a) = 5A L (a): the 
power contribution to the lattice mass Sm^ tSbt (whose perturbation theory is affected by the 
IR renormalon of the on shell mass) corresponds to the UV behaviour of the potentials while 

6 In fact this is one way to define 2<5m^ at in perturbation theory: the r-independent part of the Fourier- 
transform of the momentum space lattice potential 33 . 
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the power contribution to Ajj (whose perturbative expansion has the UV renormalon of A^ s ) 
is associated with the low energy behaviour of V 0j l. This is the same renormalon/power term 
structure as in the continuum OS/RS schemes. 

For a C r < Aqcd lattice effects become invisible and the formulae elaborated above 
apply under the replacement RS i— ► L. To illustrate this quasi- continuum limit, we eliminate 
the a~ l dependence from the expressions altogether, which is straight forward: 

E s (r) = 2m L (a) + E^+(r;a), (56) 
= 2m L (u7 1 )+V S)L {r;uJ 1 )+0(r 2 ), (57) 



where 



V a , L {r-vJ l ) = V StL (r;a)-25m^ t (a) + 25m^ t (uJ 1 ), (58) 
m L ( Z /- 1 ) = m L (a) +<5m s L tat (a) - Sm^JyJ 1 ) . (59) 

Note that the running from one scale to another is renormalon-free. For the hybrid case we 
can directly write, 

E H (r) = 2m L (z/ / - 1 ) + [V , L (r; vj 1 ) - 5A L (uf)\ + A H {u7 l ) + 0(r 2 ) , (60) 

where the combination V 0) l — 8Al replaces the V d ,rs of Eq. (JHOJ) - 

Finally, we mention that the situation r = on the lattice resembles the r <C uj l 
continuum situation. Unlike in the continuum, however, on the lattice, even at r = 0, all 
observables remain finite as a -1 provides us with a hard UV cut-off. 



4.5 Scale dependence 

As we have mentioned in the previous sections, the running of pole mass and gluelump 
energies with Uf, in the RS scheme, and with a, in the lattice scheme, is renormalon-free. 
Therefore, the functional dependence can be described by a convergent expansion in pertur- 
bation theory. Nevertheless, in order to achieve the renormalon cancellation, the same scale v 
has to be used in the perturbative expansion. This produces large logs if the scales Vf and u'j 
are widely separated and, eventually, some errors, if one works to finite order in perturbation 
theory. In the RS scheme, there exists a solution to this problem. Even though 5mRs( I/ /) suf- 
fers from the renormalon ambiguity, the difference 5mRs( z/ /) — 5^rs(^/) is renormalon-free. 
We can perform a resummation of 5m-^{vf) with any prescription to avoid the singularity 
in the Borel plane since it will cancel in the difference. We will take here the Principal Value 
(PV) prescription, which yields, 



oo 

<5 m RS (f/) = N m v f a s {v f ) c s 

s=0 



b-s 




(61) 



where 
and, 



D b (-x) = x {e- x {-x) b [T(-6) - T(-b, -x)} - cos(7r6) e - 1 T(-6)} , (62) 

POO 

T(b,x) = / dtt b - l e~\ (63) 

J X 
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denotes the incomplete T function. 

The second term in Eq. corresponds to Aj^g, once introduced in the sum of Eq. 

(jfilj) . It cancels from the combination 7 , 5m^g(^/) — 8rn^g(i/t), and we will not consider it 
any longer. The sum of Eq. (j61|) represents softer and softer singularities in the Borel plane. 
Therefore, we expect at least the difference Sm^^f) — 6m^ (i/p) to converge (although, 
obviously, we have no mathematical proof of this). Since the first three terms are known we 
can check if this actually happens. We can see that this is so with a high degree of confidence 

in Fig. nm 



-Hrs iyi) + ^rs k: 




Figure 10: — 5m^(vf) + 5rriRs(i>f) at LO (dashed line), NLO (dashed- dotted line) and NNLO 
(dotted line) according to the sum in Eq. \61\) . We take u'j- = 9.76 r^ 1 . 

We can also compare —5rn^(i , f)+5m^(i'j) with the corresponding difference, calculated 
at finite order in perturbation theory: 



-5m RS (v f ) +Sm RS ( 



2 

2 7T 



hi 



2 71 



(64) 



In 



We depict this comparison in Fig. where we take v = Vj to minimise one of the logs. We 
see how the finite order results approach the PV curve 8 , which we will use in what follows 
wherever we need the running. 

A similar behaviour holds if, instead of 5m R s, we study 5Ars- 



One may wonder if this cancellation materializes itself in practice since we only know the first three 
terms of the series. However, we checked this numerically and the results turned out to be virtually indis- 
tinguishable. 

8 For finite order computations we take a s with one, two, three etc. loop running according to the order 
in a s at which we work. If instead, we use a s with 4-loop running (the highest accuracy known until now) 
the convergence to the PV result is accelerated. 
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3 4 5 6 8 9 

Figure 11: — dm^iyf) + <5mRs(z/y) at LO (dashed line), NLO (dashed- dotted line) and NNLO 
(dotted line) in perturbation theory [see Eq. J b'4\ ) with v — i/f) versus the Principal Value 
result (solid line). We take u'j = 9.76 . 

For the lattice scheme we cannot perform an analytical resummation as higher order 
terms are unknown. On the other hand, there exist non-perturbative lattice determinations 
of the static masses [A#(a) and S(a) = A (a)] for different lattice spacings. They provide us 
with non-perturbative measurements of the running against which the finite order results can 
be tested. It is also possible to relate results in both schemes by perturbative renormalon- 
free expressions. We will investigate both, the running within the lattice scheme and the 
translation between both schemes in Sees. 15.2116.11 and 16.31 below. 



5 Phenomenological analysis of the gluelump spectrum 

We will determine the lowest gluelump energy, A#, from two different observables in two 
different schemes: from the non-perturbative difference Ejj u (r) — i? s +(r) in the continuum 
limit in the RS scheme as well as from gluelump energies A^(a) obtained at finite lattice 
spacings in a lattice scheme. Lattice and RS scheme can be translated into each other and 
we find internal consistency. We finally present results on the whole gluelump spectrum and 
compare our findings to previous literature. 

The situation discussed here is similar to the one encountered in the "binding energy" in 
static-light systems which we will address in Sec. Elbelow. These mesons very much resemble 
gluelumps, with the only difference that the source is in the fundamental representation and 
screened by a light quark rather than by a gluonic operator. 
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5.1 Determination of A^ s from the static potentials 

We intend to determine Ab from the hybrid potentials. For this purpose we will use our 
rif = lattice continuum limit data on AEjj u (r) = Eu u (r) — E^+(r) as obtained in Sec. El 
Using this difference allows us to eliminate the power divergence that appears in lattice 
simulations of the potentials (or, in the continuum OS scheme, the renormalon associated 
with the pole mass). Note that the difference has a well defined continuum limit. It is also 
interesting to see that the large distance linear term is cancelled as well. At the same time, 
Ab will still additively contribute to this combination, see Eq. (JHJ). In order to extract this 
non-perturbative constant, the perturbative difference between octet and singlet potentials 
has to be subtracted. For a reliable determination, the perturbative series has to be well 
defined and show convergence. However, this is complicated by the contribution from the 
renormalon discussed above and can only be achieved in a scheme where such renormalon 
singularities are taken into account. We have worked out the RS scheme in Sec. 14.31 which 
is well suited for this purpose. 




0.25 . 5 0.75 1 1.25 1.5 




Figure 12: Splitting between the U u and the potentials and the comparison with Eq. \6bl) 
for v = Vi [see Eq. (EJ)/ at vj = 2.5 . ^[(V^rs — K,rs)(?") + A§ s ] is plotted at tree level 
(dashed line), one-loop (dashed- dotted line), two loops (dotted line) and three loops (estimate) 
plus the leading single ultras oft log (solid line). 

We fit A B using the following equality (see Figs. Upland [121 for the quality of the fit): 

Eu u (r) - E xt (r) = Af(v f ) + V , RS (r; v f ) - V s>RS (r; Vf ) , (65) 

where the non-perturbatively obtained left hand side (Ihs) is renormalon- free but on the rhs 
the renormalon can be shifted between the two contributions, the ultrasoft matrix element 
Ab and the soft Wilson coefficient V Q — V s , at a given order of perturbation theory. This is 
why we have to specify the scheme, the RS scheme in our case, which we use to eliminate 
(or to reduce) this ambiguity. 
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Figure 13: Splitting between the U u and the potentials and the comparison with Eq. Ao'ty) 
with v = 1/r for Uf = 2.5 rjj" 1 . r [(K,RS — Vs,Rs)( r ) + A§ s ] plotted versus r at tree level 
(dashed line), one-loop (dashed- dotted line), two-loops (dotted line) and three loops (estimate) 
plus the RG expression for the ultrasoft logs (solid line). 

We fix Uf = 2.5 ro and the final result reads, 

Ag s (z// = 2.5 Tq" 1 ) = [2.25 ± 0.10(latt.) ± 0.21(th.) ± 0.08(A M g)] r x . (66) 

Note that is the only fit parameter. Also note that the above value corresponds to the 
nf = case. The errors of this determination stem from several sources (for the above fit 
we use lattice data up to distances of around 0.5 r ): 

1) "latt." denotes the statistical error of the fit: ±0.10. 

2) "th." stands for the theoretical errors. We first consider the error due to the truncation 
of the perturbative series (higher orders in perturbation theory/scale dependence). We 
obtain a first estimate by performing the perturbative expansion in a s {vj) or in a s (l/r). 
This provides us with an estimate of neglected subleading logarithms. Actually, in both 
cases one and the same number, A§ s « 2.25 r^ 1 , is obtained, which we take as our central 
value. The effects of higher orders in perturbation theory are estimated by considering the 
convergence of the determination of A^ s at each order in perturbation theory. Working with 
at s (vi), the series {2.43,2.37,2.28,2.25} is obtained. This series seems to show convergence 
for the last terms. In any case, the corrections are small. Working with a s (l/r), the series 
{2.00, 2.40, 2.31, 2.25} is obtained. This series is clearly convergent although the corrections 
are larger than when using a s (z/j) as the expansion parameter. To be conservative we will take 
the difference between the last two terms as the error made by truncating the perturbative 
series: ±0.06. There is also some source of error from the normalization constant of the 
renormalon of the singlet and octet potential. For the singlet potential (following Ref. [T5] ) 
we estimate a 10% error in Ny s , which produces a ±0.10 error. For the octet potential, the 
error is very small compared with other sources of error. Even if, conservatively, we consider 
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the general shift produced by setting 8V 0t 2 = (note that this also accounts for the error 
in perturbation theory of the octet static potential) our result only changes by ~ 0.01/0.02. 
We will neglect this error to avoid double counting. In the above analysis we have neglected 
non-perturbative effects. On general grounds they have the structure 

/y RS - V RS \ 

S np (En, - i% ) ^ r 2 F S- )+Br\ (67) 

The B r 2 term is due to r ■ O^EO type contributions in the pNRQCD Lagrangian (see Ref . jU] 
for details). The other term in Eq. ()67|) is due to r • O^ES type contributions. This produces 
a perturbative mass gap. F is the convolution of a short distance and a long distance piece, 
depending on the ratio of V^, RS — V RS over the masses of the gluelumps. For the purpose of 
estimating the uncertainty it seems reasonable to keep the lowest order in this expansion. 
This is equivalent to having a quadratic contribution, 

5 np (E Uu -E Et )^A nu ^ t r 2 . (68) 

If we introduce this term into the fit, we obtain tqKb ~ 2.30 [working with a s (l/r)] with 
A Uu _z+ ~ —OAtq 3 . We take the difference as an indication of the error due to non- 
perturbative effects. By summing linearly all the above errors we obtain ±0.21. 

3) "Ams" : this error is due to the uncertainty in = [0.602 ± 0.48] r^ 1 [SZj: ±0.08. 

We have performed the fit using lattice data within a window of inverse distances ranging 
from about z/j « 2.6 GeV down to Vj ~ 1 GeV. From the plots (see Figs. IT2l and IT3*j) one can 
actually see that the curves follow the lattice data up to values r < r . This corresponds 
to very low energies (< 500 MeV). Being conservative, we will not use data determined at 
these low energies without a better understanding of the dynamics. Nonetheless, such a fit 
would actually produce very similar numbers to the ones quoted above. This is even more 
so if a quadratic term is included. In general, introducing more lattice points reduces the 
statistical errors ("latt."). Including a quadratic term will reduce the theoretical error on 
Kb since some of the changes that occur when altering the order of perturbation theory can 
be absorbed into a variation of ^4n„-s+- However, the addition of a second fit parameter 
increases the statistical error and also the uncertainty due to Aj^g. We conclude that while 
the individual errors depend on the precise fitting details the total error remains remarkably 
stable. 

One might ask whether, in addition to A#, a reliable value of A u _ s + can be obtained. 
This however would require more lattice data at short distances as well as a more detailed 
understanding of the r 2 renormalon of the static singlet potential. 

We do not consider the S~ data in this section as we have already established in Sec. 12.41 
above that the difference with respect to the 11^ potential is proportional to r 2 to leading 
order. Hence we cannot obtain any independent new information on from these data, 
that have larger statistical errors. 

5.2 Determination of A^ s from A# 

There exists a direct determination of A.g(a) (the l + ~ or B gluelump) by Foster and 
Michael |T2] • The numerical values are displayed in Table El where we used the same r^/a 
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values as were used in this reference. It is clear from the discussion in Sec. 14.41 that these are 
perfectly sensible numbers if incorporated into a global scheme with renormalon cancellation, 
for instance, with the potentials also defined in the lattice scheme as in Sec. 14.41 In doing 
this we are able to independently determine in a different scheme. Consistency would 
require that after translating the lattice into the RS scheme the results should agree with 
each other. We will check this in this section. 





A^r /A* y r (LO) 


A* a r (NLO) 


A* y r (NNLO) 


A* s r (NNNLO*) 


A£ y r 


2.94 
5.27 
7.32 


5.33(10) 
6.99(05) 
8.36(05) 


1.59(19) 
1.97(17) 
2.21(17) 


2.82(12) 
3.20(10) 
3.55(10) 


2.37(15) 
2.88(12) 
3.25(13) 


2.41(10) 
2.89(13) 
3.16(13) 



Table 3: The inverse lattice spacing, the mass of the l +_ gluelump A^ in the lattice 
scheme, as well as its conversion to the RS scheme to different orders in perturbation theory. 
NNNLO* stands for an estimate obtained neglecting 1/N% corrections, for details see the text. 
In the last column, we state the values o/A§ s (a~ 1 ) using Eq. \6b}) and the running according 
to the PV prescription, Eq. \61)) . The errors only incorporate the statistical uncertainties as 
well as the 8 % uncertainty in Aj^r , added in quadrature, but no estimates of "theoretical" 
errors. 

The master formula that relates the lattice and the RS scheme reads (known up to 
NNLO), 

A| s (^) = A L H (a) - [5A L H (a) + 5Af(u f )} . (69) 

Both, A^ s and A^ have a power-like dependency on Uf and a -1 , respectively, but are 
renormalon-free, A^ exactly and A^- s within the precision of our estimation of the renor- 
malon contribution. This implies that the combination 5Al + SA-^s does not contain a 
renormalon either if calculated in a consistent way: 5Al(o) and 5A RS (uf) contain one and 
the very same renormalon contribution (with negative relative sign). The sum of both terms, 
expanded in terms of a s has good convergence properties (using the same normalization point 
to enforce the renormalon cancellation at each order in perturbation theory). The explicit 
expression at NNLO reads, 

5A L (a) + 5A RS (u f ) = ^a^a^u) (70) 

+ {lT17 {V2 + Vl [ ~ bl + 2(3 ° ln(z/a)]} + Vf " ^ 0i1 ) } a ^ v) + " ' ' 

where the can be found in Eqs. (|3~%j) and (J3~9"j) . &i in Eq. (pf2~|) and V 0) \ and V^i in Table El 
An estimate of the 0(af) term can be obtained from Eq. (J78j) below, under the replacements, 
Cf i — > Ca and V s>i ^ 2(V s ,i — Km)- This estimate will be subject to 0{1/Nf) corrections to 
the coefficient V3. 

In principle, Vf and a~ l need not be equal but we will take them similar to avoid large logs. 
The large numerical values of t>2 and b\ are mainly due to contributions from lattice-specific 
tadpole diagrams that arise because the breaking of Lorentz symmetry becomes particularly 
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evident at UV scales 



. This often results in badly convergent perturbative series when 



expanded in terms of ai{a). However, the convergence is vastly improved, once the series is 
re-expressed in terms of a more "physical" coupling like a s (a^ 1 ) = a£,[l— b\a 2 L /(<in)-\-- • •] (see 
e.g. Refs. [331 EE])- This is also evident from Eq. (fTOj) above as v\ ~ 3.17, (v% — biVi)/(4n) ~ 
-1.97 [and [v 3 - 2b x v 2 - (h - 2b\)vi]/(An) 2 « 14.5]. 
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Figure 14: The lowest gluelump mass A B as obtained on the lattice (diamonds), as well as 
converted into the RS scheme at NLO (squares), NNLO (pentagons) and NNNLO* (NNNLO 
estimate, circles). The error band corresponds to the result for A§ s of Eq. \6b}) . without the 
"theoretical" error, run to different scales, according to the PV prescription Eq. $61)) . The 
dashed lines, drawn to guide the eye, are explained in the text. 

We can now translate the values obtained by Foster and Michael into the RS 
scheme. The results are shown in Table 01 and are also displayed in Fig. El "NLO" and 
"NNLO" refer to translating from the lattice scheme to the RS scheme via Eq. ()70j) to 
0(a s ) and 0(a%), respectively 9 . Obviously, to leading order, A B is scheme independent. 
"NNNLO*" stand for an estimate obtained assuming that the NNNLO contribution to 5Al 
is equal to the NNNLO contribution to 5m L with the replacement of the overall factor 
Cf I— > Ca- This is correct up to 0(1/N%) effects. Finally, the conversion from the lattice to 
the RS scheme has been performed using the 4-loop running of a s at v = a -1 = Uf. This 

9 Note that the counting here differs from that used in Fig.frj in the RS scheme, where we labelled 0(a?) 
as "LO". 
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accelerates the convergence to the RS results. If, instead, we use the n-loop running of a s 
that is consistent with the order of the calculation, we still see convergence but with, in the 
NLO and NNLO cases, larger corrections. This is mainly due to the fact that within the 
present window of energies the values obtained for a s (v) from Aj^g from a one- or two-loop 
running are significantly different from those from the three-loop running (which is close 
to four- loop). The lattice prediction of that we use as an input applies to very high 
energies, such that it is important to run a s down to v > 2.5 r^ 1 as precisely as possible. 

Within present errors we can fit the data with straight lines but there will be logarithmic 
corrections and, in the gluelump data A^r , additional 0(a 2 ) = Oi^uJ 2 ) lattice artifacts. The 
figure reveals that at the lattice spacings investigated these are tiny, relative to the linear 
slope. Except for these lattice corrections the running of A^ is non-perturbatively accurate. 
Needless to say that the power dependence on a -1 is universal for all gluelumps, such that 
gluelump mass splittings have a well defined continuum limit, which is also confirmed in 
Ref. 

In lattice perturbation theory we can calculate the "running" of the gluelump data to 
0(a 2 ) [and up to O(a^) if we neglect 0(1/N 2 ) effects]. There is a renormalon ambiguity 
in the absolute value. However the slope is not affected by this. If we take the value 
A|^(7.32r _1 ) « 8.36 Tq 1 from Table EH and perform the running with NNNLO* accuracy, 
we obtain the dashed line that joins the "LO" RS(= L) points. We can see that this 
parametrization is quite close to the non-perturbatively evaluated data. Moreover, there is 
overall convergence, with higher order terms being numerically smaller in the lattice scheme. 
We will discuss this in more detail in Sec. |B] below, in the context of the static-light binding 
energy A , which has a similar perturbative expansion, up to an overall factor Cf/C^, see 
Fig. M 

In Fig. [TH we also compare the value obtained in Sec. 15.11 above [Eq. ()66j) ]. with run- 
ning according to the PV prescription Eq. (|61|). with the results obtained directly from the 
lattice determination of the gluelump mass via Eq. (|7U|) . We see clear convergence with 
alternating signs from LO (diamonds), NLO (squares), NNLO (pentagons) and NNNLO* 
(circles) towards the result calculated from the II n and E+ potentials in the previous section 
and its running (error band). Our NNNLO* estimates already agree with this error band. 
The dashed lines connecting the NLO, NNLO and NNNLO* points are the corresponding 
transformations of the curve through the LO points and just drawn to guide the eye. All 
errors displayed in Fig. HU are statistical only, plus the uncertainty on A^g. Within the 
theoretical errors of Eq. (|66|) (±0.21 r^ 1 ), in fact we already find agreement at the NNLO 
level. In Sees. 16.11 and in particular 16.31 below we will analyse the running of the binding 
energy of static-light mesons in more detail, see also Fig. 

We obtain an independent second prediction for A^ s from the gluelump data. The statis- 
tical errors are smaller in the gluelump case than those we encountered from the continuum 
potentials. In a first step we obtain the fit parameter, 

A| s (7.32r - 1 ) = [3.21 ± 0.04(latt.) ± 0.42(th.) ± 0.10(A M g)] , (71) 

from a global NNNLO* fit, 

A£(a) = A| s (^) + 
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5A L H ( 



+ 5A RS 



H 



(72) 



where we have chosen v = Vf = 7.32 r 1 . We can then convert this result into, 

Ai s (2.5r c ^ 1 ) = [2.31 ± 0.04(latt.) ± 0.33(th.)lS:li(AMff)] r o^ ( 73 ) 

using the PV running in the RS scheme. This compares well with the result from the 
potentials, Eq. (|SBJ). 

The errors displayed in Eq. (JTT|) above are due to the following sources: 

1) "latt." is the sum of the statistical error (±0.03) and the error encountered when 
varying the fit range (i.e. excluding the left-most data point): ±0.01. 

2) "th." is the sum of perturbative and non-perturbative errors. As perturbative errors 
we take the difference between NNLO and NNNLO* results (±0.20) as well as a 10 % 
uncertainty in N Va —N Vo (±0.18). To investigate possible non-perturbative effects we include 
an a 2 term into the fit. We estimate an additional ±0.04 uncertainty from this source. Adding 
these three errors linearly results in ±0.42. 

3) "Ajyfg" stands for the uncertainty due to the error of A^ro [37]: ±0.10. 

Whereas the statistical error is smaller in this determination than the one of Eq. (|66p 
and the uncertainty due to the error of is comparable in size, the systematics are less 
well under control, which is reflected in the large theoretical error. First of all, for the lattice 
gluelumps we only have the perturbative result to 0(a 2 ) with an estimate of the 0(a^) 
term while in Sec. 15.11 above we knew the O(a^) results and have an estimate of the 0(a^) 
terms. Furthermore, as the previous analysis was based on observables with a well defined 
continuum limit, we circumvented the problem of disentangling the a -1 "running" of AbTq 
from 0(Aq CD a 2 ) lattice artifacts. With gluelump data on more, and in particular finer, 
lattice spacings the latter disadvantage (which at present is however not the dominant one) 
can in principle be overcome. In conclusion, it is nice to observe perfect agreement between 
the two predictions, which enhances our confidence in the methods applied and adds further 
credibility to our error estimates. 



5.3 Higher gluelump excitations 

Now that we have fixed the energy of the lightest gluelump, we can quote absolute values for 
the remaining gluelump spectrum using the results of Foster and Michael [T2j. We display our 
predictions in Tabled where the errors correspond to the sum of the individual uncertainties, 
added linearly. The dominant uncertainty is that of A#, as the mass differences between 
the different gluelumps have been determined with very good accuracy. Needless to say that 
these results are scheme and scale dependent. The quoted numbers refer to the RS scheme 
with i>f = 2.5 Tq 1 m 1 GeV. With the information presented in this paper they can be run to 
different scales. For ease of reference we also converted these values into GeV units (using 
Tq 1 = 394 MeV). However, we note that one should add a scale uncertainty of about 10 % 
to them to account for the fact that all results have only been obtained in the quenched 
approximation. 

Note that the gluelump operators can be represented in terms of gluonic fields EH]- 
In general one and the same gluelump can be created by infinitely many different adjoint 
operators H. Within each channel we display (one of) the lowest dimensional such choice (s) 
in the table. The basic building blocks are the covariant derivative (with J PC = 1~ + , 
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H 


A$V 


A^/GeV 


1+- 


Bi 


2.25(39) 


0.87(15) 


1— 


Ei 


3.18(41) 


1.25(16) 


2~ 


B{iBj} 


3.69(42) 


1.45(17) 


2+- 


D {l E j} 


4.72(48) 


1.86(19) 


3+- 




4.72(45) 


1.86(18) 


0++ 


B 2 


5.02(46) 


1.98(18) 


4~ 


D{iDjD k Bi} 


5.41(46) 


2.13(18) 


1-+ 


(B A E)< 


5.45(51) 


2.15(20) 



Table 4: Absolute values for the gluelump masses in the continuum limit in the RS scheme 
at Vf = 2.5 Tq 1 « 1 GeV, in r units and in GeV. Note that an additional uncertainty of 
about 10 % should be added to the last column to account for the quenched approximation. 
We also display examples of creation operators H for these states. The curly braces denote 
complete symmetrization of the indices. 

dimension 1), the chromomagnetic field Bi (l +_ , dimension 2) and the chromoelectric field Ei 
(1 , dimension 2). The curl of the electric field has the quantum numbers of the magnetic 
field, such that on the lattice all states can be created by operators that are local in time. 
Furthermore, D • B and D • E can be eliminated, the first because it is identically zero, using 
the Jacobi identity, the second by applying the equations of motion. One example: the 
lowest dimensional operator that creates the 3 + ~ state is D^DjB^, where the curly braces 
denote the sum over all 10 symmetric permutations of the indices. This includes three terms 
Di J2j DjBj = such that indeed there remain only seven independent operators to create 
this seven dimensional representation. Also note that DuBj\ and DuEjx each only contain 
five independent operators, consistent with J = 2 etc.. 

It is interesting to see that the level ordering roughly corresponds to the lowest dimension 
of the creation operator, once the equations of motion are used to eliminate the E field jH]- 
This makes the E field "heavier" than a B field, increasing its dimension by one. The 3 
gluelump (two derivatives and one E which corresponds to dimension five, after substituting 
E) is not included into the table as no controlled continuum limit extrapolation was possible. 
However, its mass at fixed finite lattice spacing is in the same ball park as that of the other 
dimension five states, 4 and 1 h , in support of this naive operator counting picture. 

5.4 Comparison with previous results 

We shall relate our results to previous determinations of the gluelump masses. All these 
suffer from the problem of obtaining the global constant and, in none of these, the scheme 
was clearly defined, such that they need not yield the same results that we obtain. 

In Ref. jSH] the gluelumps were studied within a string model. One general feature of this 
approach is the excess of predicted states. This seems to be a problem of this model since it 
does not appear to be compatible with QCD, or more precisely with its realization for this 
kinematical regime: pNRQCD jHj (see also the discussion in Ref. [40 ). The prediction of 
this model, As(n/ = 0) = 1.87 GeV, is by a factor of two larger than our result. 
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In Ref. the same value for the electric and magnetic correlation length is obtained: 
A.e(n/ = 0) = A B (nj = 0) = 0.90(5)(10) GeV, from lattice simulations using the cooling 
method. The number for A# coincides with ours. However, the splitting between chromo- 
electric and -magnetic correlators is unaccounted for. From the results of Foster and Michael 
one would then assign a systematic error of the order of this splitting w 400 MeV: clearly 
a better conceptional understanding of how "cooling" removes short distance fluctuations, 
without destroying essential infrared physics would be useful. On the other hand it is com- 
forting that numbers similar to our results are obtained in this approach, that is also meant 
to subtract the perturbative contributions from the low energy matrix element. 

In Ref. a sum rule analysis of the electric and magnetic correlator was made. The 
main result was A£(n/ = 0) = (1.9 ± 0.5) GeV. It should be noted that the value of Aj^g on 
the lattice is now smaller by 5%, compared to the value used in this analysis. Taking this into 
account we find this result compatible with ours [1.25(16) GeV], within errors. Moreover, in 
this analysis, evidence for A^ > A# was reported. 

In Ref. jlH]; & n MIT bag model calculation was used to obtain the gluelump spectrum. 
No errors were assigned to this evaluation. The value of A^ is by about 500 MeV larger 
than ours and quite consistent with the sum rules evaluation. The same holds true for A#, 
however, for the higher excitations the agreement with the results of Foster and Michael is 
less convincing. 

In Ref. jllj . lattice correlation functions that are needed to calculate relativistic correc- 
tions to the static potential were used in order to check the validity of the stochastic vacuum 
model in the Gaussian approximation. Under this assumption, which was to some extent 
tested in this reference, these correlation functions could be related to gluonic field strength 
correlators and upper limits for the gluelump masses were obtained: A#(?2y = 0) < 1.64(16) 
GeV and A E (nj = 0) < 1.04(15) GeV, respectively: the ordering of the gluelumps is wrong, 
however, the upper limits quoted are in no contradiction to our results (or indeed to a 
different ordering). 

In Ref. jH], a constituent quark model was used. The results roughly agree (within a 
200 - 300 MeV error) with the splittings predicted by Michael and Foster and the hybrid 
spectrum at short distances (see Ref. jlU] for some criticism of this evaluation). For the 
lightest gluelump they obtain A# « 1.4 GeV. 

We have seen how different determinations of A# result in values ranging from less than 
one GeV up to nearly two GeV. These numbers are all scheme dependent. This may explain 
the huge differences between different results. Our result provides strong constraints on 
vacuum models. Furthermore, the RS scheme provides a unified framework to study the 
non-perturbative effects in an unambiguous and model independent way. 

6 Static-light systems 

The situation discussed above very much resembles the one that one encounters in heavy- 
light mesons in the static limit. In this case, the adjoint source is replaced by a fundamental 
source which is not screened by gluonic fields but by a light Dirac quark instead. (A light 
Higgs scalar in the fundamental representation would be an alternative possibility.) In these 
systems the binding energy A of the ~ state (which will correspond to pseudoscalar and 
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vector heavy-light mesons, once l/m^ corrections and the spin of the heavy quark are taken 
into account) plays a role similar to that of the A# discussed above. The experimental mass 
of the B meson Mb can be factorized into, 

M B = A + m b + 0(l/m b ), (74) 

where both A and m b depend on scheme and scale. In the literature (see e.g. Ref. jTH]) the 
binding energy in the lattice scheme is referred to as S(a) = A L (a), which is renormalon-free 
but has an a -1 power divergence. For the Wilson action and n/ = this <5m^ at (a) power term 
is known to 0(al) in perturbation theory [Eqs. (|37jl - (|4U|)]. Subtracting this perturbative 
result introduces renormalons. 

It is also possible to define the binding energy in an entirely non-perturbative renormalon- 
free and power-term free way, for instance by subtracting the energy of a temporal Schwinger 
line in Coulomb or Landau gauge jlSJ. In fact the same can be achieved in the case of the 
lowest gluelump mass, either by subtracting the energy of an adjoint Schwinger-line in a fixed 
gauge (see also Ref. [TH]) or by subtracting the on-shell mass of an adjoint Polyakov- Wilson 
line, encircling a compactified lattice dimension. From an EFT point of view however one 
would like to combine a non-perturbative low energy result with a perturbative calculation 
at high energies. For instance to quote a value for the b quark mass in the MS scheme, the 
UV renormalon of the binding energy is required to cancel the IR renormalon of the OS mass 
and hence a perturbative subtraction is essential: the renormalon of the expansion of the 
power divergence is the same as the one that is encountered in the conversion from the OS 
mass into the MS mass. This procedure has been implemented in the past in calculations of 
the b quark mass from lattice simulations in the static limit [T8] . 

The b quark mass has also been obtained in perturbative QCD in the RS scheme at Uf = 2 
GeV from the T(IS') system using EFTs ^Hj- Subtracting this value from the spin-averaged 
mass of the B meson yields, 

A RS (z// = 1 GeV) = [0.365 ± 0.085(th.)+°;°^(A M g)] GeV . (75) 

This number is different from the value quoted in Ref. ^H] 10 , since here we have performed 
the running to Vf — 1 GeV using the PV prescription and not included 0(1 /nib) corrections 
into the fit (these two effects partially compensate each other). Using the PV prescription 
allows us to perform the log resummation for the renormalon related terms. However, the 
result strongly depends on the value of Aj^g. 

Eq. (|75jl has been obtained from the physical T and B systems, not in the quenched 
approximation. The scale r^ 1 = 394 ± 20 MeV [2HJ EEJ EZ| is also obtained from T phe- 
nomenology. Re-expressed in terms of r we get, 

A RS (^ = 2.5^) = [0.92 ± 0.22(th.)i°;^(A M s)] . (76) 

In what follows we will extract A from lattice data of static-light mesons. After ad- 
dressing the b quark mass we will conclude with a more detailed study of the running in 
the lattice and RS schemes, using precision data from the static potential within an energy 
range, 2 < r u f = r /a < 15. 

10 Again, note that what we call RS scheme here corresponds to the RS' scheme of Ref. |15| . 
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6.1 Determination of A 

We will use Eq. ()76|) as our starting point for the n f = situation. In order to compare with 

lattice results in the quenched approximation we will employ the nj = running of A (vf) 
and keep in mind that on top of the errors stated above one might expect an additional 10 % 
quenching error. 



Ref 


crVo 


A r 


A KS r (NLO) 


A KS r (NNLO) 


A KB r (NNNLO) 


-r-RS 

A r 


[50 


2.93 


2.45( 6) 


0.79(10) 


1.34( 7) 


1.16( 8) 


0.99(1) 


M 


2.93 


2.22( 4) 


0.56( 9) 


1.11(5) 


0.93( 6) 


0.99(1) 


E3 


4.48 


2.86( 4) 


0.83( 8) 


1.37( 6) 


1.23( 6) 


1.16(3) 


E3 


5.37 


3.28( 6) 


1.03( 9) 


1.59( 7) 


1.45( 8) 


1.22(4) 


EH 


6.32 


3.44( 8) 


0.96(11) 


1.53( 9) 


1.40( 9) 


1.28(4) 


EH1 


7.36 


3.83( 8) 


1.10(11) 


1.70( 9) 


1.57(10) 


1.34(4) 


Hi 


7.36 


3.87(11) 


1.14(13) 


1.74(12) 


1.61(12) 


1.34(4) 


m 


8.49 


4.24( 8) 


1.24(11) 


1.87( 9) 


1.74(10) 


1.40(4) 


[48 


9.76 


4.49(10) 


1.20(13) 


1.85(11) 


1.72(12) 


1.45(5) 



Table 5: The inverse lattice spacing [51], the static-light binding energy A L = £ \48\ ^ 



f57Jj / in the lattice scheme, as well as its conversion to the RS scheme to different orders 
in perturbation theory. In the last column, we state the values of A^ s (a _1 ) using the PV 
running, Eq. \61\) . of the result Eq. \7b\j in the RS scheme. The errors only incorporate the 
statistical uncertainties as well as the 8 % uncertainty in Ajj^tq added in quadrature. 
The values in the last column, which have been obtained from the physical T(1S) and B 
meson masses, have additional errors inherited from Eq. \7b\ j, which, however, will only 
result in an overall upward or downward shift and will not affect their differences. 

A L (a) has been calculated on a variety of lattice spacings by different collaborations [JH1 
OH IHl EH! • The main source of uncertainty in these determinations is the extrapolation to 
zero light quark mass. We used the r /a values from the interpolation of Ref. jHIj to assign 
the scale 11 . The results are displayed in Tableland are roughly consistent with each other, 
with the exception of the coarsest lattice point ro ~ 2.93 a that corresponds to (3 = 5.7. Here 
the raw data of Ref. [201 are more accurate but the chiral extrapolation of Ref. [31] should 
be better controlled. 

We multiply the values obtained for Agr of Ref. [T3] (that are displayed in Table EJ) by 
the colour factor C//C4. At (3 = 5.7,6.0 and 6.2, respectively, we obtain the numerical val- 
ues 2.37(4), 3.11(2) and 3.72(2). The corresponding values in Table EJread 2.45(6) |2.22(4) , 
3.28(6) and 3.83(8) |3. 87(11) where for both, (3 = 5.7 and (3 = 6.2, two independent determi- 
nations exist. The qualitative agreement is remarkable: not only the perturbative expansions 
of SA and Sm sta t are dominated by terms that are proportional to the respective Casimirs 
of the gauge group representation of the static source but also the non-perturbative values 

11 These values slightly differ from those quoted in Ref. ^2 used in Table above, which cover a smaller 
window of lattice resolutions. 
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Figure 15: Perturbative running of the binding energy A in the lattice scheme, in comparison 
with lattice data, starting at the smallest available lattice spacing. The NNNLO error band 
incorporates the error due to the uncertainty in A-^g [37$ , and the statistical error. 

themselves. In fact also in the RS scheme the result Eq. (J66|) is close to the value displayed 
in Eq. (J7HJ), multiplied by C A /C f = 9/4. 

Similar to the discussion in Sec. 15.21 above, we can translate the results from the lattice 
scheme into the RS scheme. The master formula in this case is very similar to Eq. (Jfillj) and 
reads (known to NNNLO), 



A (v f ) = A (a) - <Jm£ at (a) - 8m RS {vf) 



(77) 



with 
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where, 



Bi(x) = -h + 2p i _ 1 ln(x) 



z = l,2, 



(79) 
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Figure 16: The binding energy A, obtained on the lattice and converted into the RS scheme 
at NLO (squares), NNLO (pentagons), NNNLO (circles) and the result of Eq. \81\) . run to 
different scales using the PV prescription (neglecting the "theoretical" errors). The dashed 
lines are explained in the text. 



and the coefficients V s ,i and V s ^ can be found in Tabled The coefficients V{ and 6, can be 
found in Eqs. (JSHl - (EH and Eqs. (gSJ) and (02)), respectively. 

Eqs. (J77)l and (fTSj) also relate results obtained at different lattice spacings to each other, 

AV) = A L (a) - 

To illustrate this we display the A L (a) values of TableElin Fig.HSl together with the expected 
running, starting at the finest, i.e. right-most, lattice point at LO, NLO, NNLO and NNNLO. 
The NNNLO error band contains both, the statistical error and that due to the uncertainty 
in a s (a). The running is done in each order in a self-consistent way to the given order in a s , 
according to Eq. (f78|) (without the V S) i terms). We used v = 9.76 r^ 1 and the initial value, 
« s (z/), was calculated from A^gT using the four loop running. We observe convergence 
and moreover the series is sign alternating. To NNNLO, except for the lower lying of the 
two r /a pa 2.93 data points, there is no contradiction between data and the expectation. 
However, the points of Ref. [23 have a slightly more pronounced slope such that the r Q /a ~ 
2.9,4.5 and 6.3 points (j3 = 5.7,5.9 and 6.1) lie below the curve while the r /a ~ 8.5 point 
(P = 6.3) lies somewhat above. 



^stat(a) 



^stat(a' 
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We also display the data of the table in Fig.UHl in analogy to Fig.HUbut disregard the LO 
result that is already displayed in Fig.ED The size of 5m^ at (a) — 5m~ns(a~ l ) increases linearly 
in a -1 , with logarithmic corrections: at coarse lattice spacings there might be significant 
perturbative 0(aj) and non-perturbative 0(a 2 /rl) corrections affecting the slope of this 
function while at fine lattice spacings the slope can be determined accurately but the 8m- 
difference itself becomes large. An accurate conversion between the two schemes can therefore 
neither be obtained at extremely fine nor at very coarse lattice spacings. Setting v — i/f, 
the difference between NLO and NNLO translation is minimised for 3 < r Vf < 4 while that 
between NNLO and NNNLO is minimal for 7.5 < r z/j < 9, where the widths of these bands 
are determined by our uncertainty in the value of A^gT - 

We choose to translate the lattice scheme results into the RS scheme by means of a global 
NNNLO fit to the r^/a > 5, i.e. (3 > 6.0 data, expanded in terms of ot s {y = 9.76 ), where 
we set Vf = v. The result reads, 

A RS (i/j = 9.76 r r ) = [1.70 ± 0.08(latt.) ± 0.18(th.) ± 0.04(A M g)] r^ 1 . (81) 

Note that A RS is the only fit parameter. 

The dashed curves in Fig. ^] correspond to such an NNNLO fit to the LO results, 
subsequently transformed in the same way as the data points to NLO, NNLO and NNNLO. 
The error band corresponds to the result Eq. (}8T| above, without the theoretical error, 
run to different energies, using the PV prescription, Eq. (|6"T]) : unlike the band displayed in 
Fig. E] above, this is not the result of an independent determination. We would also have 
found agreement with the result Eq. ()76|) . but only within the large theoretical errors of this 
un-quenched determination. 

At high order in the perturbative expansion and at high energies one would expect the 
slope of the non-perturbative running in the lattice scheme, translated into the RS scheme, 
to approach that of the running within the RS scheme. Discarding the four data points of 
Ref . [33] , Fig. El nicely confirms this expectation. We will investigate this running with 
higher accuracy in Sec. 16.31 below. 

The errors of the determination Eq. (|81|) above stem from the following sources: 

1) "latt." is the sum of the statistical error (±0.03) and the error encountered when 
varying the fit range a"[ n r = 4.48,5.37,6.32 (±0.05): ±0.08. 

2) "th." is the sum of perturbative and non-perturbative errors. As perturbative error 
we take the difference between NNLO and NNNLO results. Varying the fit range as above 
this difference never exceeds ±0.04. We also study the error due to the uncertainty of Ny s 
obtaining ±0.06. To investigate possible non-perturbative effects we include an a 2 term into 
the fit. We estimate an additional ±0.08 uncertainty. Adding these three errors linearly 
results in ±0.18. 

3) "Ajyfg" stands for the uncertainty in the determination of A^ro [37]: ±0.04. 

Using the running in the RS scheme (we note that the error due to the uncertainty of 
Nv s almost cancels in the running) we obtain, 

A (u f = 2.5r ) = [1.17±0.08(latt.) ±0.13(th.) ±0.09(A M g)] 1 . (82) 

from the value Eq. ()81|) . This nj = result compares reasonably well with the phenomeno- 
logical rif = 4 value of Eq. ([TH]) above and its error is of a comparable size. 
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6.2 Comment on the b quark mass 

We cannot resist the temptation to obtain a value for the RS scheme bottom quark mass, 
using Eq. (J71|) and our quenched result Eq. We obtain, 

m 6 ,RsO/ = 1 GeV) = [4849 ± 32(latt.) ± 60 (th.) ± 35 (A^)] MeV , (83) 

where we have translated Eq. (|82|) into physical units for Uf — 1 GeV and also added an 
extra theoretical error of ±30 MeV, due to l/m& corrections, combined quadratically with 
the theoretical error inherited from the lattice determination of A . From this number we 
can compute the MS scheme result, 

"\msK,ms) = t 4191 ± 29 ( latt -) ± 47 ( th ') ± ^ A Ms)] MeV > ( 84 ) 

where we have performed any running and manipulation with rif = 4 and used the PV 
prescription to run m^Rg from 1 GeV up to the bottom MS mass 12 . In this way higher 
order terms in the relation between the MS and the RS mass are minimized. If instead 
one determines ^&Ms( m fcMs) directly from its perturbative relation with m^R^lGeV) one 
obtains a somewhat larger result, but with sizeable higher order terms. Note that some of 
the theoretical errors, such as the uncertainty of Nv a , are correlated with the running of a s . 

Obviously one has to allow for quenching errors. Naively one might assume an O(10 %) 
effect on the binding energy which amounts to 50 MeV in Eq. (JH3j) . However, this might be 
an underestimate since the running of the mass with the scale in the rif = case is very 
different from that for rif = 4 and the relative effect on A < rrib, due to a different running, 
is larger than that on the quark mass. To illustrate this we also work consistently with 
nf = and obtain, 

™ b ,MsK,Ms) = [4339 ± 29(latt.) ± 49(th.) ± 9(A m )} MeV : (85) 

this differs from the value Eq. ()84|) by almost 150 MeV. Note that we have used the rif = 
value A^jg = 0.602 r^ 1 to obtain the above results. Using the rif = 5 QCD world average 
a s (M z ) = 0.118 instead (running it across the bottom flavour threshold down to 1 GeV), 
the central value of Eq. (jSIJ) would read, Ms( m fe Ms) = 4113 MeV. The difference between 
these two values may also be indicative of the typical size of the error due to quenching. 

We feel that 1 GeV might be a more natural scale to obtain an nj =4 prediction from 
the quenched model than 4 GeV and hence prefer the central value of Eq. (jH3J). After all, 
the quenched model has been adjusted to reproduce low energy QCD phenomenology and 
indeed Eqs. (I76J) and (j82j) agree with each other within errors. However, as discussed above 
and as indicated by the 150 MeV difference from using a different perturbative running, such 
predictions have to be consumed with some caution. Eq. (J84)) demonstrates the precision 
that can be achieved in lattice simulations of static-light mesons with sea quarks to NNNLO. 
Obviously, the "latt." error can systematically be reduced. Note that, with NNNLO per- 
turbative results, the dominant theoretical uncertainty (apart from the sea quark content) 
is due to 1 /nib corrections. 

12 We ignore the charm mass threshold. Since the charm quark mass is not much heavier than 1 GeV this 
is a small effect anyhow, completely paled by our dominant source of uncertainty, the rif — approximation. 
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6.3 The running of A from the static lattice potential 

To leading order, the singlet static energy E s is the sum of twice the heavy quark mass and 
the singlet potential, Eq. (J2J), while Mb is the sum of the quark mass and the binding energy 
A, Eq. (j7J). Consequently, in the OS (RS) schemes V s contains twice the leading renormalon 
(power term) of A. In QCD with sea quarks this is also evident from the large distance 
behaviour, where E s (r) will approach 2Mb- 

In the lattice scheme, the non-perturbative energy E^+ differs from E s by twice the quark 

^9 

mass, Eq. (|45|). and contains the same power term as the static-light energy A (times two). 
One can explicitly verify this in perturbation theory. In QCD with sea quarks Ey+(r) will 

approach 2A for r > r c , where r c denotes the distance associated with "string breaking" and 
is implicitly defined by, E s {r c ) = 2M B - We find the static potential HE] Ey+ to exceed 

the values of 2A of Ref. HE] at r > r c = (2.25 ± 0.15) r , a distance that is statistically 
indistinguishable from the value r c « 2.3 vq, obtained in simulations with rif — 2 light sea 
quarks E3] ■ 
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Table 6: The inverse lattice spacing, the estimate of the static-light binding energy in the 
lattice scheme, A pot = E^ + (ro)/2-\-A, Eqs. 1186]) and \8T\ j, as well as its conversion to the RS 
scheme to different orders in perturbation theory. The errors only incorporate the statistical 
uncertainties of the Ey+(r ) data, as well as the 8 % uncertainty in Ajjgr Q JffTy , added in 
quadrature. The overall error due to the uncertainty in A, which does not affect the running 
of A , is not displayed. 

The difference 2A - E^+(r ) = E^+(r c ) - Ey + (ro) is a constant, up to 0(a 2 ) lattice 

2jg 2jg 

artifacts. In what follows we will investigate the running of, 

Ap Ot (a) = i^ 9+ (r ;a) + A, (86) 

as a function of a -1 . The static lattice potential (tq)/2 can be determined more precisely 
than A L \ in terms of computer time it is cheaper to obtain with the same statistical error 
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and, since no chiral extrapolation is involved, with virtually no systematic uncertainties. 
However, we do not know the absolute normalization A. We re-analyse the lattice potentials 
of Refs. HEj; to correctly account for the propagation of the uncertainty of r into the 
combination ro-EWfro). By matching the lattice potential E^ + (ro)/2 = (2.856 ± 0.014) r^ 1 

to A = (3.844 ± 0.065) r^ 1 at (3 = 6.2 (r /a w 7.3), where we have two independent results 
for the latter quantity jlHl HH] , we obtain 



A = (0.988 ± 0.067) r^ 1 . 



(87) 



For ease of comparison with Sec. 16.11 we display the resulting A 



pot V 



in Table H as well as 



in the figures. The additional uncertainty due to the error in A should be kept in mind. 
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Figure 17: The binding energy A t; Eq. fifty), in the lattice scheme (full diamonds), in com- 



parison with A L of Sec. \6.1\ (open diamonds) . The constant A has been adjusted by requiring 
agreement between the two data sets at ro ~ 7.3 a. The uncertainty of A = (0.988±0.067) r^ 1 
is not included into the errors. NLO, NNLO and NNNLO refer to transformations of 
A pot into the RS scheme to different orders in perturbation theory. The solid line corre- 

and the central value of Eq. \81\) . 



sponds to the NNNLO expectation with A™ 



A RS ( Uf = 9.76 r^) 



0.602r^ 1 , 



1.70 r, 



-i 
o ■ 



-rL 



We display Ap Qt in Table together with conversions into the RS scheme, according to 
Eqs. ((ZZj) and (ffH|) . The data are also depicted in Fig. H3 (full diamonds), together with the 
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results from the static-light energies A (open diamonds). Except for the four data points 
of Ref. [33] at r /a ~ 2.9,4.5,6.3 and 8.5, whose slope is somewhat incompatible with the 
results from the other references as well as with perturbation theory (as we already noticed 
in Sec. 16.11 above), we find agreement between the non-perturbative running of _E E +(ro)/2 

and that of A L , down to the lowest scales. This need not be so since in principle the results 
may differ by 0(a 2 /r 2 ) lattice terms. We also compare this running with the expectation 
from the value w 0.602 r^ 1 (solid line), where we use the normalization suggested by 
the central value of Eq. (jUJ), A RS (u f = 9.76 Tq 1 ) « 1.70 r Q \ 

As can be seen there is no contradiction between the lattice data and NNNLO perturba- 
tion theory down to scales as low as 2r ( 5~ 1 and as high as 15 r^ 1 . This agreement is quantifi- 
able: a one-parameter NNNLO fit to the a -1 > hr^ 1 data (setting Ay^ = 0.602 r^ 1 ) yields 
X 2 '/N DF = 3.53/4, with the value (translated into the RS scheme for ease of comparison), 

A RS (9.76r 1 ) = (1.71 ± 0.01) Tq 1 . (88) 

j^g 

Including all available data results in x 2 /Ndf = 6.91/9 with A (9.76 ) = (1.70 ± 
0.01) Tq 1 . The errors of the above examples are purely statistical. The uncertainties in 
A and A^ as well as theoretical errors are unaccounted for. If we go to NNLO we obtain 
the x 2 /Ndf values of 16.3 (all data points), 23.0 (or 1 > 5r _1 ) and 6.7 (a -1 > 9rp Also, 

p^g 

the predicted value of A (9.76 r$ )r becomes somewhat unstable, ranging from 1.76 (all 
data points), 1.79 (a -1 > 5r ~ 1 ), up to 1.91 (ar 1 > 9r^ ): within the accuracy of the data 
it is essential to go to at least NNNLO in perturbation theory. 

In Fig. El and Table |H1 we have also displayed the results, translated into the RS scheme 
to different orders in perturbation theory. In Fig. El we focus on this comparison. This 
figure very much resembles Fig. El only that now the error bars are smaller as we discard 
the error of A, which will only affect the overall value of A but not the running with the 
scale. The dashed lines correspond to NNNLO perturbation theory in the lattice scheme 
with Aj^g = 0.602 r^ 1 , and the central value of Eq. (J88|) as normalization point. This 
running perfectly agrees with the data down to very low energies. As already observed in 
Sec. 16 .11 above, we also find nice convergence for a~ l > 3 r$ , as the order of the perturbation 
theory is increased. The error band corresponds to the PV prescription of the running 
in the RS scheme 13 with A (9.76 r^ 1 ) = (1.70 ± 0.04) r^ 1 , run to different scales, using 
A^jg = (0.602 ± 0.048) r^ 1 . Note that the errors that we display in this case are only due to 
the uncertainty in a s , with all other error sources of Eq. (|5T|) (as well as the uncertainty of 
A) ignored. 

We find excellent agreement between data and the predicted running. In fact, one can 
in principle determine a s from the logarithmic corrections to the a~ l running of the binding 
energy: in dedicated lattice simulations of the short distance static potential tremendous 
statistical accuracy can be achieved and tiny lattice spacings are accessible [23]. Even using 
our static singlet potentials [321 ESI that are ^ ess accurate than those of this recent reference, 
a two-parameter NNNLO fit to the a' 1 > 5r^ data yields, Aj^ = (0.590 ± 0.036) r^ 1 

13 At Vf 3> v = 9.76 Tq 1 we find some differences between the NNNLO running in the lattice scheme 
(dashed black line) and the PV prediction (error band), due to large logs in the difference Eq. I)8()[l. where 
we have not attempted a log resummation. 
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Figure 18: The binding energy Ap Qt; Eq. A86\) . translated into the RS scheme at NLO 
(squares), NNLO (pentagons) and NNNLO (circles). We have neglected an overall error 
in the vertical scale of ±0.067 r^ 1 , due to the uncertainty of A, that does not affect the run- 
ning. The dashed lines correspond to the NNNLO running in the lattice scheme with A-^as = 
0.602 r$ l , where we used the fit result Eq. A (9.76r _1 ) ps 1.71 , as normalization. 

^ j^g 

The error band corresponds to the prediction Eq. i fffT)) . A (9.76 ) = (1.70±0.04) r$ , and 
includes the uncertainty due to Aj^g = (0.602 ± 0.048) r^ 1 (but no other errors). 

and A RS (9.76 r^ 1 ) = (1.73 ± 0.04) r^ 1 with x 2 /Ndf = 3.35/3. Including the whole energy 
range, down to a -1 ~ 2rg , results in, = (0.627 ± 0.026) r^ 1 and A (9.76 r^ 1 ) = 
(1.68 ± 0.02) ro , still with very acceptable x 2 /Ndf = 6.08/8. The results for A are in 
perfect agreement with those obtained in Eqs. (JHTj) and (jHHj) above. Moreover, the fits are 
consistent with the value of Ref. [3*7j . A^ = (0.602 ± 0.048) r^ 1 , within statistical errors 
smaller than the uncertainty of this reference value. 

In conclusion, we have demonstrated that the running of the binding energy in the lattice 
scheme can be reproduced with incredible accuracy in NNNLO perturbation theory, in terms 
of a s . This accuracy is possible since, unlike in the case of the binding energy itself, there 
is no leading renormalon contribution to its running. Down to energies of about 1 GeV 
we do not see any sign of a break-down of perturbation theory or evidence of significant 
non-perturbative contributions to the running. We have also confirmed that the theoretical 
errors estimated in Eqs. (JSTj) and (J%2J) are indeed conservative. 
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7 Gluinonium and other related issues 



We already mentioned that gluelumps are interesting in the context of bound states including 
heavy adjoint particles, such as gluinos of SUSY models (even if it is quite likely that they 
will decay before any kind of hadronization takes place). In this case, to leading order in 
HGET (Heavy Gluino Effective Theory), the gluino mass can be obtained from the relation, 

M d = A° s + m 5 ,os = Al» + m §iKS (iy), (89) 

in a scheme of choice that can then be converted into the mass in say the MS scheme rrig(rng), 
analogously to the discussion of Sec. HI above. We will limit most of the discussion below to 
the RS and OS schemes but translation into lattice schemes is straight forward. 

Mq denotes the mass of the lightest (spin-averaged) glueballino. Note that in this context 
the gluelump energy A# plays the same role as the binding energy A did for heavy-light 
mesons. We have A^f{vf) = A^ s — 5Ars(z/^) and hence, 

™>gfls(yf) = m ~g,OS + SA RS (p f ) : (90) 

$ m g,RS — — m t ne glueballino case corresponds to the Sm^s of heavy-light mesons. We 
can also write down the above equations in the lattice scheme in which case, using the same 
conventions as in other parts of this paper, Srrig t L = SA L . 

In addition to glueballinos one can imagine gluino-gluino bound states: gluinonia, T. 
Their dynamics is dictated by the following Lagrangian, 

£ P NRQCD,r = / d 3 r d 3 RTr |s r + (id - V At .) S r + Or,! 1 " («A> - V A , ) O r ,i + •••}, (91) 

at leading order in 1/rrig and in the multipole expansion. This is analogous to Eq. ((TJ), 
replacing the static sources in the fundamental by static sources in the adjoint representation. 
This means that there will be further multiplets in Eq. (|91j) that we will not consider in this 
paper. 

The singlet potential between two adjoint sources V A , s {r) has been calculated in pertur- 
bation theory to 0(a 2 ) and the corresponding energy E^(r;a) was determined in lattice 
simulations (see e.g. Ref. [21] )• Up to lattice artifacts oc a 2 we can write, 

E A (r) = 2m hL {a) + E L A {r ] a) (92) 
= 2m g:L (a) + V A , s>L (r;a)+0(r 2 ) (93) 
= 2m~ g , RS (is f ) + V As , RS (r;is f ) + 0(r 2 ), (94) 

where the normalization of E^(r) with respect to E A (r) can be obtained from the gluinonium 
spectrum. Obviously, 

limE A (r) = 2M 6 , (95) 

while for the bottomonium energy in QCD with sea quarks one obtains (up to l/m^ correc- 
tions and neglecting radial and gluonic excitations), 

lim EJr) = lim E H (r) = 2M B . (96) 
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In combining Eq. (}9"4*|) with Eqs. (ffitj) and (JHUj) we obtain the important equality, 

E A {r) + 2[E B {r)-E s {r)\ = V A;S , RS (r;u f ) + 2[V , E s(r;iy f )-V S! ns(r-,u f )] + 2M d + 0(r 2 ) , (97) 

where we have used the fact that Mq = m^Rs + A§ s and E B G {Eu u , E^-}. The effect of 
Smjis cancels from the combination E B — E s and 5Ars from E A + 2Eh- Since the glueballino 
mass is a physical observable this implies that, up to 0(r 2 ) corrections, the combination 
Va,s(?~) + 2[V (r) — K(r)] is scale independent and free of renormalon and power contribu- 
tions: the UV renormalon of V Q is cancelled by the UV renormalon of V s while the leading 
IR renormalon of V , which we studied in this paper, is cancelled by one half of the UV 
renormalon of V AjS . in fact, to 0(a 2 ) this combination explicitly vanishes and the 0(a^) 
term is suppressed by a colour factor 1/N 2 . 

In the above equation E B (r) corresponds to the U u or S~ hybrid levels. For a general 
Eh{t) we would have to replace the Mq on the rhs by the mass of the excited glueballino in 
the respective channel. At r — > oo the rhs of Eq. (}9~Tj) will approach 2Mq, see Eqs. (J95J) and 

dm. 

We wish to compare our expectation with lattice data. This can either be done after 
an extrapolation of these to the continuum limit or at finite lattice spacing in the lattice 
scheme. Re-expressing Eq. (19 7|) in terms of the energy levels as determined in the lattice 
scheme [E A (r) = E A (r; a) + 2m~ g ^{a) etc.], and using the conventions of Sec. 14.41 above, this 
amounts to, 

E L A {r; a) + 2[E^( r] a) - E^(r; a) - A£(a)] (98) 

= V A>s , L (r; a) + 2[V , L (r; a) - 5A L (a) - V s , L (r; a)} + 0(r 2 ) + • • • 
= 0(r 2 ) + 0{al/Nl) + 0(A 2 QCD a 2 ) + 0(a 2 /r 2 ) . 

Both Ihs and rhs of the above equation are explicitly free of a -1 power terms (and of leading 
renormalons) . In fact the rhs vanishes in perturbation theory, to at least 0(a^/N 2 ). As 
indicated in the equation, in general there will be non-perturbative 0(a 2 /r 2 ) as well as 
0(AQ CD a 2 ) lattice artifacts, in addition to the 0(r 2 ) corrections from higher terms in the 
multipole expansion. 

The above combination is extremely interesting as for small r there should only be a 
quadratic but no linear term. At r > 2 r the adjoint string will break and the Ihs of the 
equation will approach zero like 1/r. In the intermediate region 0.5 ro < r < 2r one 
would expect two non-perturbative contributions, a linear term from the slope of E A {r), 
with an effective string tension [23], a e f[ = [3.09 ± O.lOjr^ 2 , as well as a 1/r term that 
dominantly originates from the difference between static hybrid and singlet potentials and 
whose coefficient will approach 2tt as r — ► oo, in an effective string model expectation. In 
fact for r m r one would expect this 1/r term still to dominate over the linear term. 

We wish to compare this expectation to numerical data. Unfortunately, on isotropic 
lattices where we know the gluelump mass in the lattice scheme we did not compute the 
adjoint potential while on our anisotropic data sets all potentials, singlet, adjoint and hybrid 
are available but the gluelump mass is unknown. In Fig. we display the combination 
E\(r) + 2[E^ u (r) - E^+(r)] as a function of r = r[l + 0(a 2 /r 2 )], Eq. (fTUjl . at our finest 

anisotropic lattice spacing, a a ~ 0.16ro ~ 4a r which, within errors, is compatible with the 
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Figure 19: The combination E\(r) + 2[E-^(r) — E^ + (r)}, Eq. \98\). as a function of r = 

r[l + 0(a 2 /r 2 )], Eq. MU\) . together with various fits, on an anisotropic lattice with resolution 
a a m 4a T m 0.16 r . 



continuum limit, see Sec. Eland Ref. [21]. Note that there is an additional 1 % overall error 
on ordinate and abscissa due to the conversion from lattice units into units of ro that we do 
not display. 

From Eq. (J55J) we would expect the combination shown to approach the gluelump energy 
in the lattice scheme, A^(a), as r — > 0. We see that the approach towards this limit 
is remarkably flat. In fact, excluding the r > 0.9 r data, which are clearly in the non- 
perturbative regime anyway, we are unable to resolve deviations of the data from a constant. 
Note that the units on the ordinate, 0.2 r^ 1 « 80 MeV, are quite small. A linear plus 
quadratic fit, 

E L A {r) + 2[E£>) - E^(r)} = 2A L B + cr 2 , (99) 

to r < 0.5 r data yields, 

2A| = (15.51 ±0.10) r^, c = (0.07 ± 0.70) r 3 . (100) 
A purely phenomenological fit to the same functional form for all distances results in, 

2A^ = (15.45 ±0.06) r \ c = (0.38 ± 0.07) r^ 3 , (101) 
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while in a physically completely unmotivated funnel parametrisation, 2A^ + e/r + kr, we 
obtain, 

2Ajj = (14.95 ± 0.20) r^ 1 , e = (0.08 ± 0.04), k = (0.84 ± 0.18) r 2 : (102) 

the r dependence is so weak that on the 1 % error level of the lattice data we are unable to 
discriminate between different parametrizations. However, we can determine the gluelump 
mass rather precisely, A^(a CT , a T ) = (7.75±0.05±0.07) r^ 1 , where the second error reflects the 
uncertainty in the lattice determination of r§ja T . In fact the same can be done for the a CT ~ 
0.23 r and a a ~ 0.33 tq data sets. The respective results read, = (6.71 ± 0.04 ± 0.09) r^ 1 
and A|? = (5.75 ± 0.10 ± 0.05) r^ 1 , respectively. The data are in agreement with a linear 
slope in a" 1 but, unfortunately, at present we only know the NLO perturbation theory for 
the anisotropic case. After subtracting twice these gluelump energies, we find scaling of the 
coarse lattice data with the results depicted in the figure, within error bars of comparable 
size. 

In particular from the fit to the funnel type parametrisation we see that the data leave 
little room for perturbation theory style short- distance Coulomb terms. This is in agreement 
with our expectation. However, miraculously there is also no evidence for a quadratic term 
in the r < 0.9 ro data and in fact we can set the limits 0.46 > cr^ > —0.18 for such 
a contribution, from the r < 0.7 tq data. We believe that the smallness of this term is 
accidental as had we replaced the n n by the £~ potential, it would certainly be present, 
see Sec. 12.41 One can however speculate that there might be a cancellation of r 2 effects and 
that U u does not receive a significant r 2 contribution in the multipole expansion. This issue 
should be addressed in future theoretical and numerical studies with enhanced accuracy. 

The observed slope at larger distances (~ 0.84 r^ 2 ) is much smaller than that of the 
adjoint potential in this region (^ 3.09 r^ 2 ), in agreement with our expectation that the 1/r 
contribution from the difference 2(Ejj u — E^.-) cannot be neglected. 

There is no evidence of a linear non-standard short- distance term for r < 0.9 r , at least 
not of the size expected in various models [53]. A possible explanation of the absence of 
such a term from our calculation of a quantity that vanishes to low orders in perturbation 
theory would be that a s (q) itself receives 0(l/g 2 ) corrections (see Refs. [SB])- We remark 
to this end that a s is not a physical observable. In the MS scheme it is perturbatively 
defined. The difference between a s and any non-perturbative generalisation of this coupling, 
that would allow inclusion of such singularities, will necessarily not be universal but depend 
on the prescription used. However, we are investigating a physical observable here that is 
scheme independent. 

Combinations of different potentials that lead to renormalon and low order perturbation 
theory cancellations are certainly an arena worthwhile to study for a determination of higher 
order terms in the multipole expansion and for testing the validity of the standard operator 
product expansion picture. As we shall detail below many such combinations exist. 

There are also hybrid excitations in the adjoint channel. The perturbation theory in this 
case is richer than for potentials between fundamental sources as8®8 = l©8©8©10© 
10* © 27: in addition to singlet and octet, we have another octet, two decuplet fields and 
a 27-let which have to be included into Eq. (j9*T|) . Consequently, adjoint hybrid potentials 
cannot only have the octet perturbative expansion but some will correspond to decuplets and 
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others to 27-lets. Note that the decuplet representation is not self-adjoint but has vanishing 
triality as it should be. 

The renormalon of the octet potential between adjoint sources is the same as in the 
fundamental case but the decuplet and 27-plet adjoint potentials contain new renormalons, 
which are related to those of the singlet potentials between colour charges in these respective 
representations. This exactly resembles the situation discussed above where the adjoint sin- 
glet potential contains the same renormalon as the fundamental octet potential. In fact one 
can define an infinite tower of states with different renormalons following this construction, 
a theoretically interesting enterprise but not likely to be of much direct phenomenological 
use. 

The inclusion of the octet states of Eq. (pTTj) is necessary for any consistent pNRQCD 
calculation of gluino pair production near threshold at NLO [H3|. At NNLO the decuplet 
and 27-plet fields will also play a role. In fact such contributions, depending on the mass 
of the gluino (and on its existence), might be of bigger importance than in the case of ti 
production because there are more of them. This is an exciting and very clean-cut situation 
since v 2 and r _1 are bigger by a factor ~ Ca/C/ than for quarkonia, such that all "soft" 
physics is clearly and extremely safely within the perturbative domain. 

Let us finally mention that A#, the binding energy of the lightest glueballino, determines 
the size of the splitting between the adjoint singlet potential and the lowest adjoint hybrid 
potential at short distances, the latter of which, unfortunately, has never been determined 
in lattice simulations. This is very different from the case of fundamental sources where 
binding energies of heavy-light systems, A, are much smaller than the gaps between ground 
state and hybrid excitations. In "hadrinos", that contain stable adjoint sources, gluonic 
excitations would hence play a very prominent role and simple constituent-gluino models 
might fail terribly. Unfortunately, in nature we do not encounter such particles. It would 
however be most entertaining to confirm this expectation in lattice simulations. 

8 Conclusions 

We report compelling evidence that for distances around 1 GeV -1 the gluonic excitations of 
the static potential are in the short distance regime. 

We are able to obtain a value for the lowest lying mass A# of the bilocal gluonic correlation 
functions with well controlled uncertainties, by fitting to the difference between the II U and 
£+ potentials. The RS scheme result reads, 

A| s (n/ = 0) = [2.25 ± 0.10(latt.) ± 0.21(th.) ± 0.08(A M g)] r Q l , (103) 

for vj = 2.5 Tq 1 w 1 GeV. Translated into physical units this reads, 

A| s (l GeV) = [0.887 ± 0.039(latt.) ± 0.083(th.) ± 0.032(A M g)] GeV . (104) 

Note that one should also add an extra error of order 10 % due to quenching to these 
numbers. With the information presented in this paper v$ can be run to different scales (see 
Fig. H4*)l . We also obtain values for the masses of other gluelumps, listed in Table HJ as well 
as for the non-perturbative slope ^4 nii _ E - = 0.92to;52 r^ 3 °f the quadratic difference between 
the lowest two hybrid potentials. 
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In order to state sensible numbers for A#, the scheme for the renormalon cancellation 
has to be specified. Otherwise, very different numbers can be obtained, as we can see from 
a comparison of the result in the lattice and the RS schemes. One can translate from one 
scheme into the other in a renormalon-free way, order by order in perturbation theory and 
check whether both results are consistent with each other. We have been able to confirm 
this. If we use the gluelump results from Foster and Michael at finite lattice spacing 
we obtain, A| s (2.5r c T 1 ) = [2.31 ± 0.04(latt.) ± 0.33(th.)ioi9(A5s)l ^ which is perfectly 
compatible with the result Eq. (jl(J3|) above, albeit with slightly larger errors. 

We also investigate the binding energy of heavy-light mesons in the static limit and to 
NNNLO in the conversion. We arrive at similar conclusions. For the binding energy we 
obtain the rif = value, 

A (u f = 2.5 r ) = [1.17±0.08(latt.) ±0.13(th.) ±0.09(A M g)] 1 , (105) 

which is in good agreement with the phenomenological value, obtained from the experimental 

T(1S) and B meson systems |TS], A RS = [0.92 ± 0.22(th.)±o:li (^ms)] r o 1 - 

We have demonstrated the internal consistency of our approach. Lattice predictions 
for A^ and A at different lattice spacings have been studied. We have shown that the 
perturbative series, Eq. (jS0J), relating A^(a) and A L (a) with A^(a') and A L (a'), respectively, 
in the lattice scheme is free of renormalon singularities and has nice convergence properties, 
as indicated by the consistency with the non-perturbatively obtained values. In particular 
this means that from the knowledge of A^ and A at a given lattice spacing values at 
different lattice spacings can accurately be predicted. We have studied the conversion of 
lattice predictions for A^ and A into the RS scheme. This conversion is also dictated 
by a perturbative series which is free of renormalons. We have verified that the values in 
the lattice scheme indeed approach the results in the RS scheme with a convergent pattern 
and, remarkably, the ^/-scale dependence predicted by the RS scheme is reproduced, within 
errors. We remark that for the z//-scale running it is possible to obtain a resummed non- 
perturbative expression in which the renormalon is cancelled and at the same time the log 
resummation is performed. 

We stress that the RS scheme used here is designed to smoothly converge to (MS-style 
dimensional regularization) perturbation theory at low orders in a s ; after all, the renormalon 
effect only sets in asymptotically at large orders in perturbation theory. Different values for 
As and for A can be obtained in other schemes but only at the inconvenience of having 
large corrections to "standard" perturbation theory at low orders. In this sense we consider 
our approach "natural" ; the RS scheme incorporates salient features of both, dimensional 
and lattice regularization. The approach readily benefits from results computed in the MS 
scheme, the scheme in which perturbative quantities are usually known to the highest order. 
On the other hand, by subtracting renormalons we encounter explicit power divergencies, 
which is exactly what one obtains with a hard lattice cut-off too. 

Our model independent non-perturbative predictions can directly be incorporated into 
perturbative calculations, within effective field theories, or exploited in the context of QCD 
vacuum models or calculations based on non-local condensates. Obvious phenomenological 
applications in the context of EFTs are pNRQCD in the kinematic domain mv 2 < Aqcd < 
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mv, translating glueballino masses into RS or MS gluino masses within HGET (Heavy Gluino 
Effective Theory), or converting heavy-light meson masses into quark masses within HQET. 

We observe that A B ps (Ca/Cj)A « mc/2, where mc denotes the mass of the lightest 
glueball. The first similarity is not necessarily surprising since there are technical parallels 
between A#, which corresponds to the binding energy of an adjoint source, and A, the 
energy of a fundamental source. We do not intend to advocate a constituent gluon picture. 
Nevertheless, it may seem reasonable that the binding energy of the glue to an adjoint 
source has about half the size of the energy of an entirely gluonic state. It should however 
be noted that the latter is an unambiguously defined state in the physical spectrum while 
for the binding energy A# we necessarily encounter the scheme and scale dependence that 
we discussed. 

We have also investigated the scenario of gluinonia and other excitations in non-funda- 
mental channels. While gluinos might not exist in nature and certainly do not form light 
bound states, such that phenomenological applications are limited, from a theoretical and 
conceptual point of view the existence of this part of the spectrum is very interesting. The in- 
clusion of such potentials allows one to identify many combinations in which renormalons and 
other un-wanted contributions vanish, opening up a window to the study of non-perturbative 
short distance physics. 
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